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ABSTRACT 

We present foUowup optical g' , r', and i' imaging and spectroscopy of serendipitous X-ray sources 
detected in 6 archival Chandra images included in the Chandra Multiwavelength Project (ChaMP). 
Of the 486 X-ray sources detected between 3 x 10^^^ and 2 x 10~^^ (with a median flux of 3 x 
10~^^) erg cm~^ s~^, we find optical counterparts for 377 (78%), or 335 (68%) counting only unique 
counterparts. We present spectroscopic classifications for 125 objects, representing 75% of sources 
with r* < 21 optical counterparts (63% to r* — 22). Of all classified objects, 63 (50%) are broad 
line AGN, which tend to be blue in (g* — r*) colors. X-ray information efhciently segregates these 
quasars from stars, which otherwise strongly overlap in these SDSS colors until z > 3.5. We identify 
28 sources (22%) as galaxies that show narrow emission lines, while 22 (18%) are absorption line 
galaxies. Eight galaxies lacking broad line emission have X-ray luminosities that require they host an 
AGN (logLx > 43). Half of these have hard X-ray emission suggesting that high gas columns obscure 
both the X-ray continuum and the broad emission line regions. We find objects in our sample that 
show signs of X-ray or optical absorption, or both, but with no strong evidence that these properties 
are coupled. ChaMP's deep X-ray and optical imaging enable multiband selection of small and/or 
high-redshift groups and clusters. In these 6 fields we have discovered 3 new clusters of galaxies, two 
with z > 0.4, and one with photometric evidence for a similar redshift. 
Subject headings: galaxies: active - surveys - X-rays: galaxies - quasars: general 



1. INTRODUCTION 
1.1. X-ray Surveys and the Cosmic X-ray Background 

' X-ray surveys provide fvmdamental advances in our knowledge of the X-ray universe and indeed the universe as a 
whole (e.g., the Einstein Medium Sensitivity Survey (Stockc ct al. 1991); the Cambridge-Ca mbr idge ROSAT Serendip- 
itous Su rvev llBovle et al . 1997); ROSAT International X-ray/Optical Survey ijPaee et alJl996j) . the ASCA Large Sky 

, Survey ijAkivama et al.ll200Q,l . ROSAT (0.1-2.4keV) and more recently Chandra (0.3-8 keV) have resolved 80-90% of 
the Cosmic X-ray background (CXRB) into discrete sources (Hasinger et al. 1998; Rosati et al. 2002; Moretti et 
al. 2003) most of which are the unobscured AGN familiar from optical and soft X-ray surveys. However, the high 

' energy spectrum of the CXRB is much harder (F ^ lA-^) than that of known AGN (F ^ 1.9). Population synthesis 
models that satisfy both the CXRB spectrum and the observed X-ray number counts vs. flux relation (logA-log^; 
Comastri et al. 1995; Hasinger et al. 1998; Tozzi et al. 2001) favor absorbed AGN as a dominant component of the 
CXR B. In these models, most of the accretion luminosity in the universe is from obscured sources (fHasinger et alJ 
ian fc Iwasawal[T999|) . which appear to have hard X-ray spectra because circumnuclear gas absorbs low en- 
ergy X-rays. X-ray spectral analyses and opti cal fo l low-u p of fain t hard X-ray sources d etected by Chandra and XMM 
have confirmed this interpretation generally (|Brandt et aLif200lt iAlexander et alJl2001fl but also presented some sur- 
prises. Type 2 acti ve galact ic nuclei are mostly found at z < 1, and the required sp ace density of such objects must be 
much greater (Rosa tiet al.i r2002: Ale xander et al.lf200lHHornschemeier et aLll2001l) than those of standard unabsorbed 
broad-line AGN (BLAGN). For the standard AGN unification model to survive, wherein optical Type 1 and Type 2 
classifications represent different viewing angles on identical objects (Antonucci 1993), it must begin to encompass 
the population of X-ray absorbed AGN and its evolution in number density, covering fraction and/or optical depth of 

absorbers. 

Absorption may be increasing with luminos ity or redshift (Elvis et aLl ll998l : iR.eeves et al.lll99^ . and could be as- 
sociated with early circumnuclear starbursts ('Guainazz^£^^L|j200QO . The intriguing suggestions of a z ~ 1 peak in 
X-ray-selected (X-S) galaxies (Bargcr et al. 2002; Toz zrer^Ll I^^ ^ has been sugge sted as eviden ce for an epoch of 
enhanced activity related to the assembly of massive galaxies fFra nceschini. Braito. & Fa dda 2002). At hi gh redshift 
(z > 4), a sig nificant dropoff in the co- nroving space d ensity of quasars seen in optical (e.g., Schmidt 1995; Warren et all 
Il99l lOsmedlT98^ and radio surveys ijShaver et alJll99^ hints at either the detection of the onset of accretion onto 

Electronic address: |pgreen@cfa.harvard.edu| 

^ r is the photon number index of an assumed power-law continuum such that Ne{E) = Neq ■ In terms of a spectral index a from 
fu = fiyo'^°'y define F = (1 — a). 
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supermassive black holes, or a missed high-redshift population, possibly due to intrinsic absorption. Based on pre- 
liminary evide nce for constant sp ace densities of X-ray selected quasars beyond a redshift of 2, deep (ROSAT) soft 
X-ray surveys ijMivaii et al.ll20dol) have been used to support the latter interpretation. Unfortunately, the sample size 
is small with only 8 quasars beyond a redshift of 3. At these early epochs, higher rates of galaxy interactions and 
mergers are expected to have triggered nuclear activity in galaxies (e.g. Blain et al. 1999, Ostcrbrock 1993). As 
gas-rich protogalaxies grow by merging, they may induce growth in the central black holes. Preliminary models tie 
together AGN evolution wi t h hier archical growth of clustering (e.g., Cole & Kaiser 1989) and galaxy formation models. 
iWilman. Fabian fc Nu1s"ct? ('200(/) use the assumption that seed blocks for clustering each contain black holes of mass 
1.6 X 1 0^ Mg^. D uring these early high accret ion phases, quasars may be self-clo aking, naturally copious producers of 
metals l|Kuraszkiewicz et al.ll2002f) and dust l)Elvis. Marengo, fc Ka rovska 2 00^ . A truncation of growth in baryonic 
mass is per haps achieved by th e wind-driven gas from the AGN, yielding today's observed black hole/bulge mass 
correlation ijMerrit fc Ferrares3 l2'0011. 

At more recent epochs, the bright optically-selected (0-S) quasar population appears to fall dramatically between 
z ~ 2.5 and the present, on a timescale of about 2Gyr, quite rapid compared to the fading of galaxy formation in 
hierarchical clustering models. Strong luminosity evolution dominates the luminosity function (LP) during this epoch 
across radio, optical, and X-rays, possil ply associated with depletion of gas reservoirs by accretion of smaller group 
companions ijCavaliere fc Vittorinill2000|) . 

1.2. The need for the Chandra Multiwavelength Project (ChaMP) 

To understand the formation and evolution of accretion onto supermassive black holes, and its links to galaxy 
formation, an accurate census of accreting objects is required. Effectively, we must account for the accretion radiation 
energy density across cosmic time consistent with the mass density of supermassive black holes in the Universe (e.g., 
Yu & Tremaine 2002; Elvis, Risaliti fc Zamorani 2002; Mainieri et al. 2002). Galaxies without a dominant active 
nucleus also produce X-rays from by-products of star formation such as low- and high-mass X-ray binaries (LMXBs 
and HMXBs), supernova remnants, and hot diffuse gas. The distinction between star-formation and AGN-induced 
emission is a subject of long debate that relates intimately to the above discussion. Progress on the composition of the 
CXRB, and the accretion and star formation history of the universe all require a wide area, multiwavelength survey 
with greater sensitivity and completeness than previously achieved, and with reduced bias against absorbed objects. 

Both soft X-ray and optical surveys suffer strong selection effects due to intrinsic obscuration and the intervening 
Lya forest. Consistent with the hypothesis that the CXRB is dominated by absorbed AGN, a substantial popula tion of 
reddened AGN is being found i n IR (Gregg et al. 2002; Cutri 2001 ^asci. Dr inkwatcr. fc Webster 1999: Wilkes et alJ 
I2OOI and radio ijBecker et al.i ri997: Masci 1999; Ves ter et al.ll200lF surveys. While radio surveys are least affected 
by intrinsic absorption, radio-loud objects constitute a minority of AGN. Optical classification rests on secondary 
properties (e.g. emission lines from ionized plasma). By contrast, the primary signature of accretion is emission from 
close (10-100 gravitational radii) to the supermassive black hole, which is why only sensitive hard X-ray selection 
samples all known varieties of AGN. Furthermore, for Chandra and XMM, sensitivity to emission up to 10 keV 
observed fraine) can reveal hidden po pulations of active galactic nuclei (AGN) including heavily obscured quasars 
Norman et 111120021 IStern et al.ll2001() . Absorption (~ 1022-23 cm-2) has 10 - 3000x less impact on detection by 
Chandra (0.3-8 keV) than on ROSAT (0.3-2.4 keV). High-z objects can be detected through an even larger intrinsic 
absorbing column of gas and dust because the observed-frame X-ray bandpass corresponds to higher energy, and hence 
more penetrating X-rays at the source. 2 Therefore, optical and X-ray surveys will complement each other, providing 
a fair census of mass accretion onto black holes at high redshift. 

With its low background and tight PSF (<,1" FWHM on axis) Chandra is probing the high energy universe with 
unparalleled sensitivity and resolution. The Chandra Multiwavelength Project (ChaMP) is a wide-area study of a 
large sample of serendipitous X-ray sources detected in Chandra archival images. The ChaMP has selected more 
than 100 Chandra images, in which we will detect several thousand serendipitous (non-target) extragalactic sources . 
A first X-ray catalog of these sources and early results are published in accompanying papers l)Kim et al.ll2003afR) . 
The ChaMP AGN sample includes faint nearby AGN, absorbed quasars to 2; ~ 3 and unobscured AGN to z ~ 5, 
populating new regions of L — z space, and providing a broad census of the population that will reap statistically 
significant samples for rare source types. The two 0.1deg2 Chandra IMscc Deep Fields probe far too little volume 
at low redshift where many of the absorbed AGN probably reside (Barger et al. 2001; Fiore et al. 2000) and contain 
few bright examples to easily characterize these sources. The ChaMP will effectively bridge the gap between flux 
limits achieved with the Chandra deep field (CDF) observations and those of past surveys. The flux regime where the 
ChaMP has greatest s ky area (10~^^ erg c m~2 s""'^) coincides with the break betw een power-law slop es fitting the hard 
(2-10 keV) logA^logS l|R,osati et al.l l2002^. Furthermore, as shown in Figure 1 o f|K im "^^11 12003^ a gap m coverage 
remains between deep ASCA and Chandra surveys near the break. The wide-area and medium depth of Chandra fields 
included in the ChaMP will span these fluxes with good statistics. 

A number of new high galactic latitude serendipitous surveys are under way both with the XMM-Newton X-ray ob- 
servatory and Chandra. The XMM-Newton Survey Science Center (Watson 200ll) targets several flux ra nges, including 
a med ium sensitivity survey (/90.5-4.5keV)> 2 x 10~^'*erg cm~2 s~^) with first results described bv iBarcons et all 
l)2002j) . The XMM EPIC instrument covers a large (^ 15' radius) field of view and has sensitivity out to higher 
energies than Chandra. However, given the higher background and larger PSF of XMM (5-7" FWHM), the minimum 

^ The observed-frame, effective absorbing column is A''^^^ ~ Nn/{1 + z)^'^ (Wilman & Fabian 1999). 
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detectable flux for an on-axis point source with a typical AGN spectrum is similar. For cross-correlation to faint 
sources in other wavebends, the tight Chandra PSF is a great advantage. 

While the Chandra Multiwavelength Project (ChaMP) emphasizes study of the nature and evolution of active 
galaxies, as sketched above, it will also provide to the scientific community well-defined wide-area survey products 
suitable for detailed multiwavelength studies of quasars, stars, star formation, galaxies, clusters of galaxies, and 
clustering. 

We have chosen to analyze a subsample of 6 ChaMP fields for this paper as an example of what can be expected 
from the full ChaMP sample. After describing in 32 the overall ChaMP field selection process, we brief ly describe 
the pr ocessing and analysis of the X-ray images in which is detailed in our companion X-ray paper ijKim et alJ 
l2003a^. We then outline our ongoing optical imaging campaign ([Q, X-ray/optical source matching ([0, and optical 
spectroscopy (50. We describe overall sample results in with separate emphasis on AGN, galaxies, and clusters of 
galaxies. Prospects and plans for ChaMP science are sketched in fJHl We discuss selected individual objects in each of 
the 6 fields in Appendix 1. Throughout this paper, we assume Ho=70 km s~^ Mpc~^, f^A = 0.7, and JIm ~ 0.3. 

2. FIELD SELECTION AND X-RAY IMAGE PROCESSING 

Our field selection, and X-ray image processing strategies are detailed in the companion X-ray paper of iKim et alJ 
((20023)1 Briefly, fields are restricted to high Galactic latitudes {\b\ > 20°). We omit a handful of pointings from other 
serendipitous surveys like the Chandra Deep Fields, and exclude fields with large bright optical or X-ray sources. This 
selection yields 146 ACIS exposures in Cycles 1 & 2, with Chandra exposure times ranging from 2 to 190 ksec, and 
covering about 14deg^. The sky placement of these fields is shown in an equatorial Aitoff projection in Figure ^ The 

full field list is public, available at 

|\protecl||ittp://hea- www.harvard.edu/CHAMP/ 




Fig. 1. — Location of all 137 ChaMP high latitude (|fe| > 20dcg) fields is shown in an Aitoff equatorial projection, with RA and Dec in 
decimal degrees. Filled circles represent ACIS-I at the aimpoint, and open circles ACIS-S. The average area for ChaMP-selected ObsIDs 
corresponds to 5.8 CCDs turned on, each of which is 8.3' on a side. Point size crudely indicates Chandra exposure times, varying from 2 
to 190 ksec. 



4 



Green et al. 



Table 1. Chandra Fields 



ObsID 




i_J^lJ*^0 lAL vl. 


ACTS CCDs'' 


RA 


DEC 


TIT Datr 








(kscc) 




J2000'^ 




(102"cm~2) 


536 


MS1137.5+6625 


114.6 


012567 


11:40:23.3 


+66:08:42.0 


1999 Sep 30 


1.18 


541 


V1416+4446 


29.8 


01237 


14:16:28.8 


+44:46:40.8 


1999 Dec 02 


1.24 


624 


LP944-20 


40.9 


23678 


03:39:34.7 


-35:25:50.0 


1999 Dec 15 


1.44 


861 


Q2345+007 


65.0 


23678 


23:48:19.6 


+00:57:21.1 


2000 Jun 27 


3.81 


914 


CLJ0542-4100 


48.7 


01257 


05:42:50.2 


-41:00:06.9 


2000 Jul 26 


3.59 


928 


MS2137-2340 


29.1 


23678 


21:40:12.7 


-23:39:27.0 


1999 Nov 18 


3.57 



''Effective screened exposure time for the on-axis chip. BI chips 5 and 7 generally have lower screened exposures since they are more 
susceptible to solar flares. 

''The ACTS CCD chips used in the observation, with the aimpoint chip in italics. The ACIS has 10 chips total, of which CCDs 5 and 7 
are back-illuminated. Due to telemetry limitations, at most 6 chips can be read out. 
•^Nominal target position, not including any Chandra pointing offsets. 

The fields studied in this pilot paper include 6 early (Cycle 1) fields listed in Tabled for which we have already 
reduced and calibrated both archival X-ray images and photometric NO AO 4meter/Mosaic imaging. These 6 fields 
represent about 12% of ChaMP fields of similar exposure time or greater. Four of the 6 fields have clusters as their 
intended Chandra targets. We do not include the intended (PI) targets in our analysis, their presence may bias the 
sample of source types analyzed here. The large number of fields in the full ChaMP and the variety of PI target 
types allow importa nt tests for possible sample biases. W hile an excess of point sources in cluster fields has been 
claimed anecdotally ijCappi et al. II2OOII iMartini et al.ll2002D. no significant difference is found using the larger sample 
of Chandra images (29 with clusters and 33 without) in lKim et alTl|2003bD . 

3. X-RAY IMAGE PROCESSING AND ANALYSIS 

We have used Chandra archival data reprocessed (in Apri l 2001) at CXC'^ . The foundation of the ChaMP 
rests on our own Chandra image analysis pipeline XPIPE ijKim et a n I20n3aj^ . built mainly with CIAO tools 
(fVprotect attp://asc. harvard.edu/ciao), and designed to unifor mly and ca r efully screen and analyze the X-ray data. 
Details of XPIPE can be found in our X-ray companion paper. iKim et all l|2003ar . Briefiy, XPIPE screens high par- 
ticle background periods, cosmic rays and hot pixels from each ACIS CCD exposure. Wavelet transform detection is 
performed (using wavdetect; Freeman et al. 2002) at 7 scales from 0.5 to 32" thus accommodating the Chandra PSF 
variation as a function of off-axis angle from the mirror axis. We select a significance threshold parameter of 10~^ per 
pixel, which corresponds to about one spurious pixel detection per CCD. We have confirmed this rate by an extensive 
simulation showing that the number of spurious sources detected is always less than this rate. An additional flux 
criterion of signal to noise (■^) greater than 2 from our X-ray photometry further decreases the number of spurious 
sources in the sample. We avoid spurious sources near detector edges by generating an exposure map for each CCD 
with the appropriate aspect histogram, and then applying a minimum exposure threshold of 10% to all pixels included 
for analysis. 

The absolute celestial positions of sources detected with wavdetect are accurate to about 1", which is the specified 
absolute position accuracy for the Chandra observatory."* To check for and minimize any systematic offset in the 
absolute source positions returned by wavdetect, the positions of X-ray sources are correlated to optical sources in the 
field. The average position offset from X-ray sources to matched optical sources is then applied to each X-ray source 
position. This procedure is described in more detail in 30 

We visually inspect all the X-ray images to validate sources and flag them for contamination (by hot pixels, 
readout streaks, bad bias values) or excessive uncertainty due to contamination of source or background regions. 
Throughout this paper and for ChaMP statistics, the target sources are flagged and excluded from analysis in all 
fields. Source catalogs and FIT S images are being made publicly available, with details at the ChaMP web site, 
|\protecl||ittp: / /hea- www-harva rd . edu / CHAMP / 

At the" X -ray source positions derived from wavdetect, we perform aperture photometry. We use a circular source 
aperture whose radius encompasses 95% of deposited energy, determined at 1.5 keV, but restrict this radius to the 
range 3" < rap < 40". We use a background annulus extending from 2 to 5 times rap, excluding flux within rap of 
any other detected sources in the annulus. We perform photometry in 3 energy bands: soft (5; 0.3-2.5 keV), hard {H; 
2.5-8 keV), and broad {B; 0.3-8 keV). We calculate net count rates using the effective exposure and vignetting factors 
both at the source and background regions. 

Only sources with ^ > 2 in at least one of the bands B, S, or H are included for analysis here.^ Flux is calculated 
from the net count rate in the B, or if necessary from another {H or S) band to satisfy our criterion. We assume a 

3 CXCDS versions R4CU5UPD14.1 or later, along with ACIS calibration data from the Chandra CALDB 2.0b. 
* See \protect ittp://asc. harvard.edu/cal/ASPECT/ for a detailed discussion of Chandra aspect solutions. 

^ is calculated from cr = cr'^^^ + ^^'^^kg ^^^re atot is the error in the total number of counts in the source aperture, and cri,i^g is the 
error on counts in the source-excised background aperture. R is the ratio of areas of source and background regions. The jr criterion is a 
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powerlaw slopes from F = 1.0 to 2.5 causes (for Ng""^ = 3 X lO^Ocm-^) an increase in the derived /x (0.5-2 keV) of 
~ 12% for FI chips, and 6% for BI chips. 

The input bandpass energy range is passed to PIMMS (Mukai 1993; v3.2d, which includes Chandra Cycle 4 effective 
area curves) along with the adjusted count rate and spectral model parameters to calculate a de-absorbed flux fx in 
the 0.5-2 keV band, since this bandpass is most useful for comparison to numerous published results from ROSAT, 
XMM, and Chandra. The error in fx is computed by applying the same counts-to-flux conversion factor to the error 
in counts, increased by a further 15% to qualitatively account for errors induced by the spectral model assumptions. 
These results are listed on-line along with other quantities described below. 

For every source, we wish to calculate an effective on-axis ACIS-I hardness ratio 



that may be fairly compared among sources to study their X-ray properties even in the low counts regime. Many 
(especially weak) sources have derived H or S counts that are formally < 2(T, or even negative. Naive inclusion of 
these yields nonsensical HR values outside the physical range from -1 to 1. Therefore, we first replace such counts by 
a 2(7 upper limit. This yields for HR either a 'detection' (when both H and S counts are > 2a), an upper limit (when 
H < 2a h), a lower limit (when S < 2a s), or very rarely an indeterminate value (both). 

Mirror reflectance changes as a function of energy and of the angle 9 between the source and Chandra^ s aimpoint, 
causing energy-dependent vignetting. A different energy dependence arises from the detecting CCD - the quantum 
efficiency of backside- (BI) vs. frontside-illuminated (FI) chips. Therefore, for an identical source spectrum, different 
HR values result, depending on 9 and chip. To correct for this, we perform a second adjustment of H and S counts, 
normalizing them to refer uniformly to an on-axis source falling on an FI chip. To do this, we have modified the 
PIMMS code to include vignetting. Where the integration of source model is convolved with the instrument effective 
area (EA), we incorporate an extra multiplicative factor - the percentage of the on-axis EA as a function of photon 
energy and 9. We interpolate between analytic fits (courtesy Diab Jerius) to EA as a function of 9 from to 30', 
performed for 4 energies (0.49, 1, 1.48, 2.02, 2.99, and 6.4 keV). Similar curves are shown in the CXC Proposer's Guide 
Rev4.0 fCXC v4.0 2001). 

To derive the corrected hardness ratio HRq, we first convert counts to de-absorbed flux in each band using our 
standard source model, the source 9 and chip. Using the same source model, we then convert the derived flux back 
counts at 9 ~ on an FI chip. The resulting corrections in HRq do not depend strongly on the assumed F; the 
difference in the mean HRq derived assuming F — 0.5 and 3.0, is only 0.02 (3%) out to 15' off-axis. 

We have performed a series of simulations to derive the expected HRq values for high sources across a grid of 

F and Nh values, as shown in Figure [3 A comparison to the median corrected value of HRq = —0.5 for all sources 
detected in both H and S, and assuming no absorption implies F ~ 1.4. 
Finally, if at least one of H and S are significant, then we calculate errors for HRq from the H and S counts (or 



Optical imaging is crucial to provide optical fluxes, preliminary source classification, and accurate centroiding for 
spectroscopic followup. Optical centroids supersede X-ray centroids for extended or multiple objects, for cluster 
galaxies, or for Chandra sources with large off-axis angles and hence less accurate X-ray positions. Multiband colors 
can provide source classifications and photometric redshifts for the majority of fainter ChaMP sources (r* > 22) for 
which high quality spectra are more difficult to obtain. 



For every ChaMP field, we scale our optical exposure times to the Chandra X-ray exposure times, to provide a 
uniform sensitivity to X-ray/optical flux ratios. Our color-limited survey thus minimizes optical telescope usage, while 
accessing for every fleld a similar fraction of every object type. To convert X-ray counts to flux, we flrst assume 80% 
of the exposure time remains after cleaning high particle background periods.^ The low background and small PSF 
of Chandra allows high confidence detection of point sources with as few as 4-5 photons, but to reduce uncertainties 
on fx and thereby improve object class discrimination, we adopt 10 counts as the ChaMP minimum detection 
limit. We convert the detectable number of counts for each Chandra field to an X-ray flux limit, and then scale the 
corresponding optical magnitu de limit to include > 90% of ROSAT AGN even at the X-ray flux limit, based on 1448 
ROSAT-detected quasars from lYuan et al.| ljl998jl . The typical resulting magnitude limit as a function of Chandra 

flux, not a detection criterion, since detection is based on probability in wavdetect. Sources with = 2 can be easily detected in many 
cases, since the background rates are typically 1 broadband (0.5-10 kcV) count per 40 pixels in 100 ksec ICXC v4.0 ll20()lll . 
® Since X-ray photons are tagged with arrival times, high background periods can be excised. 
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Fig. 2. — Simulated on-axis FI chip Hardness Ratios for high sources. Top. The family of curves show the resulting HR 
for different (z = 0) absorbing columns as a function of power-law Photon Index F. Bottom. The family of curves show the 
resulting HR for different values of power-law Photon Index F as a function of absorbing column. 



exposure time T can be expressed as y = 14T7 -f 2.5 logT. This criterion (logy^ < 0.9) should include /ar(7er fractions 
of other source types with brighter relative optical emission like stars and most galaxies. Sources that are relatively 
weak optically compared to ROSAT quasars will have a lower completeness and a brighter effective limiting magnitude 
in our survey (e.g., sources suffering from substantial optical extinction). Details of our optical imaging for these 6 
fields are presented in Tabled] The resulting desired limits range from V mags (note that V ^ {g' + r')/2) of 21 to 
26 with a median of 23.5 mag for all ChaMP fields. However, to achieve reasonable limits on the amount of required 
optical imaging time, we limit our maximum desired optic al magnitude limit to ~ 25 in any field. 

NOAO 4-meter imaging with the Mosaic CCD cameras (|Muller et alJll998f) is key to the ChaMP, since it provides 
adequate depth, spatial resolution (~ 0.26"/pixel), and a large field of view (36' x 36') over the full Chandra FoV. For 
shallower northern fields, for secondary calibration of deep imaging from non-photometric 4 meter (4m) nights, and 
for imaging bright objects within deep fields, the ChaMP also uses the FLWO 1.2m with the 4shooter CCD on Mt 
Hopkins (0.32"/pixel, 22' field). In this paper, we report only deep 4m fields observed in photometric conditions. 

We dither the 4m/Mosaic pointings to allow for better sky flat construction, averaging over defects, cosmic ray 
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Table 2. Optical Imaging 





Hiyr)- V ) 


lelescope 




r liter 


Uitners 


rljxposure 


Airmass 


r VV xllVi 


1^ 

niTO 
















(totcil sec) 




( I 


Limit 


Limit 


536 


0.0131 


KPNO 4m 


12 Jun 2000 


g' 


2 


1200 


1.39 


2.2 


22.875 


23.9 










r' 


2 


1200 


1.43 


2.0 


22.875 


24.2 










i' 


2 


1200 


1.48 


1.9 


22.625 


23.8 


541 


0.008 


KPNO 4m 


12 Jun 2000 


g' 


2 


1000 


1.33 


1.9 


22.875 


23.9 










r' 


1 


500 


1.40 


2.0 


22.125 


23.7 










i' 


1 


500 


1.50 


1.7 


22.125 


23.5 


624 


0.014 


CTIO 4m 


28 Sep 2000 


g' 


3 


810 


1.05 


1.2 


24.875 


25.9 










r' 


3 


630 


1.01 


1.1 


24.125 


25.3 










i' 


3 


720 


1.00 


1.1 


23.375 


24.4 


861 


0.0246 


CTIO 4m 


28 Sep 2000 


g' 


3 


1260 


1.23 


1.4 


24.625 


25.8 










r' 


3 


1080 


1.17 


1.4 


24.375 


25.3 










i' 


3 


1170 


1.17 


1.1 


23.375 


24.3 


914 


0.0383 


CTIO 4m 


28 Sep 2000 


g' 


3 


990 


1.04 


1.2 


25.125 


26.1 










r' 


3 


810 


1.06 


1.2 


24.375 


25.4 










i' 


3 


900 


1.09 


1.2 


23.375 


24.2 


928 


0.051 


CTIO 4m 


28 Sep 2000 


g' 


3 


900 


1.01 


1.6 


24.375 


25.6 










r' 


3 


720 


1.02 


1.3 


24.125 


25.2 










i' 


3 


810 


1.05 


1.0 


23.375 


24.3 



^FWHM of point sources in final stacked images. 

''Turnover magnitude limit at ~ 90% completeness, using 0.25 mag bins before extinction correction, as described in the text. 
Magnitude limit for a ~ 5fT detection. 



removal, and gap coverage. While the optical data reduction is conceptually straightforward, it is in practice compli- 
cated by the size of the data set (285 MB per image), and by the need to interactively manage bad pixel flagging. 
After photometric analysis, our final reduced images are stored in (short) integer format, which reduces their size by 
a factor of two, with < 1% effect on measured magnitudes. 

4.2. Filter Choice 

The ChaMP uses Sloan Digital Sky Survey (SDSS) g' , r' , and i' filters ijFukugita et al .''1996^. whose steep trans- 
mission cutoffs improve object classification and photometric redshift determination ((Gonzalez fc Maccarond l2003 
iRichards et 111120(11 . While is clearly our primary discriminant of AGN from stars, at least 2 optical filters are 

required to also constrain stellar spectral class. St ill, with just 2 filters, high-z, reddened or obscured AGN may share 
the stellar locus with M dwarfs l)Lenz et al.lll998() . so we chose 3 filters to separate such quasars from stars because 
these types of AGN represent a prime science goal of the ChaMP. Furthermore, it has already bee n demonstrated 
that AGN of redshift 3.5 < z < 5 can be constrained to Az ~ 0.2 with SDSS g\r\ and i' colors l|Richards et alJ 
I2OOI . Since g' and r' also bracket the 4000A break in 2 > 0.8 cluster ellipticals, this aids distance estimation and 
cluster membership evaluation. The bright A5577 night sky line falls between the g' and r' filters, allowing fainter 
limits to be achieved in shorter integration times. Finally, someday the need for NOAO northern imaging foUowup 
will be relieved by the SDSS itself for the shallower^ ChaMP fields. The huge database of SDSS magnitudes and 
colors help with deeper fields and core science as well by offering (1) bootstrap calibration of deep ChaMP imaging 
obser ved in non-photometric co nditions (2) cali bration and testing of accurate photometric redshifts (Richards et al] 
l2?)0l and source classification ((Newberd 11999(1 (3) detailed comparisons between SDSS optical and ChaMP X-ray 
sample selection. 

4.3. Optical Image Reduction 

The reduction of the raw optical CCD images is primarily done using standard techniques in IRAF.^ The Mosaic 
cameras utilize the multi-extension FITS format, allowing the ei ght individual CCD images to be treated as a single 
image for most reduction purposes. The IRAF package mscred (" Valded [20021) is used to properly handle the multi- 
extension FITS images and to simplify the data reduction. Although an overview is pre sented here, a detailed reduction 
walkthrough is available ^ for the NOAO Deep Wide Field Survey (NDWFS) data ((Jannuzi fc Devlll999() using the 
same instruments. 

Standard calibration images are taken at the telescope (bias, dome flats, etc.) excluding twilight sky flats. We 
construct a master super-sky flat by combining multiple object frames, thereby rejecting real objects in the frame 
and leaving us with a high image of "blank" sky. These improve the flat-field correction provided by conventional 

The SDSS achieves ~ 10 to g'<22mag at best IRichards et al.ll2001^ . 
* IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci ence Foundation. 
^ jhttp: / /www.noao.edu/noao/noaodeep/ReductionOpt/frames.html| 
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dome flats, mostly by accounting for differences in each filter between the night-sky and our dome lamp color. Before 
the super-sky flat can be made, images from the KPNO 4m also require subtraction of a pupil ghost caused by light 
back-scattered from the telescope optics, which affects primarily the inner four CCDs. Because the pupil ghost changes 
with the amount of light entering the telescope, its strength varies with the filter used, and with ambient light from 
the moon or bright stars near the field of view. For this reason we create a template pupil ghost scperately for each 
filter, taking care to exclude images with extremely bright saturated objects near the center of the field of view. Once 
a template pupil has been generated, it must be scaled and subtracted from individual object frames in each filter, a 
process that is now largely automated. 

4.4. Astrometry and Photometry 

We first derive 0.3" RMS) astrometry from our optical imaging to enable accurate X-ray source matches, and 
optimal positioning of spectroscopic fibers or slits. The (J2000) position of the optical sources are referenced to the 
Guide Star Catalog II. Final positions used by the ChaMP for all matched sources are the average of positions 
measured from all opt ical filters in which an o bject is detected. 

We use SExtractor l|Bertin fc Arnout^ll996D to detect sources, and measure their positions and magnitudes. A first 
run of SExtractor is used to measure the mode F of the FWHM of stellar objects in each image. For this step, our 
filters first remove objects that are unreliable (any with processing error flag > 3 - cosmic rays, spurious blended, or 
contaminated sources, or sources near a chip edge), sources that are too faint (typically < 2 x lO"* net counts) or too 
bright (typically > 2 x 10^ counts), and extended sources (SExtractor stellarity < median). For program fields, 
we typically use a minimum detection threshold (SExtractor detect_thresh) of 0.8ct above background per pixel, and 
require a minimum grouping (detectjninarea) of r^/4 pixels at this threshold for a detection. Since we are cross- 
correlating to X-ray sources - rare on the sky compared to optical sources - we are somewhat more concerned with 
completeness than with spurious source rejection. From detailed comparisons of dithered and deep stacked images of 
the same fields, we find that these parameter settings forge a good compromise between detecting most real sources, 
but not too many spurious sources. 

A second run of SExtractor measures 3 instrumental magnitudes: one from a circular aperture with diameter FWHM 
(for highest see Howell et al. 1989), another with 4 times the FWHM (which we find encompasses 95% flux for 
point sources), and flnally SExtractor 's mag_ auto. The latter is a elliptical aperture magnitude fala lKronlll980l) that 
is robust to seeing variations, as described in Ngnin o et alJ (|1999fl . Following the convention of the early data release 



of the SDSS quasar catalog l)Schneider et al. 12002 ^. we present the optical photometry here as g* , r* and i* since 



the SDSS photometry system is not yet flnalized and the NOAO fllters are n ot a perfect match to the SDSS fllters. 
Standard stars calibrated for the SDSS are available from fSmith et all l)2002|) but are generally too few per fleld to 
practically allow for a complete photometric calibrati on, especially gi ven the long readout times for the Mosaic. To 
facilitate the use of the many stand ard stars listed bv |Landoli ^ | l)1992D . we transform them to the SDSS photometric 
system using iFukueita et al.l l)1996j) to derive our nightly photometric solutions. In the standard star flclds, we flnd 
standard stars automatically by positional matching of accurate coordinates for Landolt standards (Hcndcn 2002) to 
SExtractor source positions, using a high (^^ 7a) SExtractor detection threshold to provide unambiguous detections. 

We merge instrumental photometry results from each filter into a multicolor file by positional matching using a 
search radius of 1". If multiple matches are found within this radius, the closest match is retained. If no matching 
detection is found in some filter, the RMS of the background for that filter at that object's position is recorded for 
later conversion to a limiting magnitude. In these 6 fields, our merged catalog of Mosaic photometry contains about 
343,000 objects with flags < 4 indicating good photometry in at least one band. 

We then perform photometric calibration of the standard stars for each night via a 2-step multilinear method l)Hardiel 
^12) J adapted from an implementation written in IDL by James Higdon. For a given night, for the instrumental 
magnitude v in each filter, we solve the linear equation 

vq - V = {ko + kiX) + e(co) -I- z 
where vq and cq are the standard star's cataloged magnitude and color, respectively, e the color coefficient, z a zeropoint 
correction, and fcg and ki are extinction coefficients. We iteratively solve for these 4 coefficients, first freezing fco and 
fci and solving for e and z^ then freezing e and z and solving for fco and fci, until the solution converges (typically 
3 or 4 iterations). We then delete measurements with residuals exceeding 3cr, and repeat the solution until no such 
residuals remain. 

Similarly, for each color c we iteratively solve the linear equation 

Co = e(c + fco + kiX) + z 

for the four coefficients e, z, fco and fci, deleting measurements with residuals exceeding 3cr as before. We tested the final 
results from our procedure against that oflHarris. Fitzgerald. & Reed (1981) as implemented in the IRAF photcal 
procedure, and find results identical to within the errors. Our calibrations produce e, z, fcO, fci coefficients for each of 
(7*, r*, i*, {g* — r*)^ and (r* — i*). Since these calibrations for these NOAO fllters have not been published previously, 
we list them in Table O as a reference. The final transformations are applied to the standard stars to derive RMS 
residuals for each night for g* , r*, i* , {g* - r*), and (r* — i*). These are fisted in TableOlfor each night. The nightly 
transformations applied to our program objects are 

c* = e(c + fco + kiX) + z 

The Guide Star Catalogue-II is a joint project of the Space Telescope Science Institute and the Osservatorio Astronomico di Torino. 



The ChaMP: Optical FoUowup 9 



Table 3. Nightly Photometric Calibration Coefficients 



Fit Variable 


e 


z 


fco 




N Stars 


Fit RMS 




(mag) 










(mag) 






KPNO 4m 12 June 2000 








0.060 


25.314 


0.284 


-0.217 


52 


0.057 




-0.055 


25.424 


0.207 


-0.155 


49 


0.068 




-0.185 


25.194 


0.191 


-0.117 


54 


0.250 




1.021 


-0.087 


0.044 


-0.036 


54 


0.111 


r*-i* 


0.978 


0.362 


0.064 


-0.050 


47 


0.031 



CTIO 4m 28 September 2000 





0.038 


25.695 


0.257 


-0.210 


43 


0.027 


j.*b 


-0.020 


25.794 


0.145 


-0.118 


43 


0.033 


j*b 


0.083 


25.358 


0.089 


-0.072 


41 


0.041 


3*-r* 


1.011 


-0.087 


0.110 


-0.090 


41 


0.037 


r*-i* 


0.903 


0.410 


0.078 


-0.064 


43 


0.027 



''Uses (g*-r*) colors in fit. 
''Uses (r*— i*) colors in fit. 



V* = v^(kQ^ kiX) + e[cj) + z 

where c* and v* are the final calibrated color and magnitude, respectively. 

For objects not detected in any optical filter, we assign a magnitude limit 77150- corresponding to a flux of Scr detection 
where a is the background RMS at that position. Starting fi:om the standard CCD equation, we find 

, fDTa\ 
niDa = -2.5 \ogy — — — J + 2; 

Here, D is the desired f . The RMS noise at the position of an object on the CCD is denoted as cr; F is again the 
FWHM in pixels, T is the total exposure time, and z indicates that the photometric calibration is applied, for which we 
assume the mean color for objects in the field. This limiting point source magnitude is listed for each field in Table [3 
for the median background RMS in the field. Objects at these magnitudes are not detected with high completeness. 
From our comparisons of individual dithered and deep stacked images, we find that the magnitude where the number 
counts peak in a differential (0.25 mag bin) number counts histogram corresponds approximately to 90% completeness 
in the magnitude range 20 - 25. This turnover magnitude 777to is typically about 1 mag brighter than the 5cr limiting 
magnitude. Both magnitude limits arc listed in Table [3 

We perform several sanity checks on the final photometry. First, we plot color-color diagrams against a crude mean 
stellar locus from the SDSS Early Data Release (EDR). The locus we derive by eye (from objects with < O.lmag error 
among 30, 000 high latitude objects) i n the SDSS EDR corresponds to one line primarily following the SDSS halo/thick 
disk locus iChenll200HfYasuda et alJl20?nl) 

(r* - i*) - 0.4029(g* - r*) - 0.007 

and another following the thin disk locus 

(r* - i*) - 14.286(g* - r*) - 18.5 

An example of these color loci is shown in Figure EH to be discussed in more detail in § 17.31 By studying a variety of 
SDSS fields at 7 different galactic latitudes (6 = 30, 35, 45, 60, -45, -55, -65) we find variations of about 0.2mag in 
(r' — i') occur in the thin disk locus. 

We also compare our magnitudes in g* and r* to magnitudes derived from the GSC2.2. We have derived the 
transformation between SDSS EDR and GSC2.2 magnitudes by direct comparison of the two surveys, using only stars 
between 17th and 19th magnitude. After 3tT clipping, the typical residuals are 0.2 mag, with zeropoints differing by 
- 1%. 

g* = 1.032 J + 0.1255( J - F) - 0.851 

r* = 0.982 F -f 0.1518( J - F) + 0.451 

i* = 0.994F- 0.1235(J-F) + 0.281 

We note that most optical point sources brighter than 18th mag suffer from saturation effects in our NOAO 4m 
imaging. While these objects constitute a small fraction 6%) of our optically identified X-ray sources, they are a 
larger fraction (^ 10%) of our spectroscopically-classified objects. We expect a full release of the SDSS will substantially 
alleviate this problem, as will our ongoing imaging study of these fields with the SAO FLW01.2m and 4-shooter CCD 
on Mt Hopkins. 

Since most of our X-ray sources are extragalactic, before calculations of optical luminositie s, we correct the optical 
magnitude measurements for Galactic extinction using the maps of .Schlegel. Finkbcincr. fc Davi^ llT99^ . We assume 
an Rv — 3.1 absorbing medium, and absorption in the SDSS bands as given bv fSchneider et all l|2002|) 
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5. X-RAY TO OPTICAL SOURCE MATCHING 

For the ChaMP, iKim et all i2003al) have carried out extensive simulations of point sources generated using the 
SAOSAC raytrace program (\protect ittp://hea-www. harvard.edu/MST/) and detected using CIAO/wavdetect 
(Freeman et al. 2002). For weak sources of ~ 20 counts between 8-10' off-axis from the aim point, the reported 
X-ray centroid position is correct within 1.8", corresponding to a la confidence contour. 

The source naming convention of the ChaMP as officially registered with the lAU is given with a prefix CXOMP 
(Chandra X-ray Observatory Multiwavelength Project) and affixed with the truncated J2000 position of the X-ray 
source (CXOMP Jhhmmss.siddmmss) after a mean field offset correction is applied, derived from the positional 
matching of optical and X-ray sources in each field. 

5.1. Matching criteria 
5.1.1. Automated Matching 

We obtain candidate optical identifications of X-ray sources in each field using a 2-stage automated position matching 
process between optical and X-ray source positions. In the first stage, the subset of X-ray sources within 5' of the 
Chandra aimpoint and with X-ray B (0.3-8 keV) detections > 20 counts are matched to optical sources from our CCD 
imaging, using position matching with a radial offset matching tolerance, in arcsec, iJmatch — ^/logio(-B)i where 9 is 
the off-axis angle in arcminutes, and B is the X-ray source counts, but i?match is limited to a minimum of 3" for all 
sources. This matching tolerance encompasses the wavdetect centroid positional uncertainty variatio n with off-axis 
angle and source counts, using the refined source position estimator described in IKim et al.l l|2003aV The optical 
position errors are neglegible in comparision. In those cases where there are multiple optical matches, the closest 
match is used. The median X-ray to optical position translational offset of these matched sources is then applied to all 
the X-ray sources for the field, to bring the optical and X-ray source lists onto the same coordinate frame. The X-ray 
position corrections are typically less than 1", which is consistent with the specified absolute position accuracy^ ^ for 
the Chandra observatory and wavdetect source positions. Figure |31 shows optical matches overlaid directly on the 
ACIS-I image of ObsID 536. 

In the second stage of position matching, optical matching is performed for all X-ray sources with X-ray B detections 
> 5 counts, using the same matching tolerance i?match as above. Multiple optical matches for an X-ray source are 
allowed in this stage, and are passed to the visual inspection process described next. Extended X-ray source centroids 
may have larger position errors than are predicted by our point source formula. Furthermore, in some cases no 
corresponding optical cluster galaxy may exist close to the cluster X-ray centroid. If no optical counterpart is found 
in the automated matching routine for these or other cases, visual inspection provides a second chance to identify 
counterparts. 

5.1.2. Visual Inspection 

After our automated 0/X matching between detected sources from XPIPE and SExtractor source positions has been 
performed, we examine the images by eye. An example of the plots used for visual inspection is shown in Fig 01 In 
many cases, no obvious optical counterpart can be detected, while in other cases there are several (source confusion). 
Even with typical astrometric accuracy of an arcsec, at these faint optical magnitudes and X-ray fluxes we frequently 
find several optical counterpart candidates. Beyond r* ~ 20, galaxies dominate the optic al number counts, and at by 
24th magnitude, cumulative source counts are near 10^ deg~^ l)Kiimmel fc Wagner 1120011) . At these faint magnitudes, 
the probability of finding an unrelated optical source by chance within any 3" circle becomes ~ 20%. The small PSF 
and accurate centroiding of Chandra become crucial at these faint magnitudes, as does careful source inspection. 

5.2. Fraction with Optical Counterparts 

During visual inspection of the optical field of each X-ray source, we assign a match confidence to each plausible 
counterpart. To consider a sample as uncontaminated as possible by incorrect matches, we have removed from further 
optical analysis below 45 X-ray sources with poor or ambiguous counterpart assignments. 

The resulting fraction of X-ray sources with confidently matched optical counterparts is shown in Tabled We have 
not yet achieved the desired optical magnitude limits (set to match optical counterparts for 90% of ROSAT AGN) in 
all these deep Chandra fields, but the imaging has served to facilitate spectroscopy to r* ~ 21, and to provide fiuxes 
and colors to still deeper limits. Lower achieved match fractions also reflect the inclusion by Chandra of new 4^ 

Jopt 

populations atypical for ROSAT such as heavily absorbed AGN. 

Excluding the 10% of sources with more than 200 X-ray counts, the mean/median X-ray flux for matched sources 
(3.5/2.4 in units of 10~^^ erg cm~^ s~^) exceeds the flux for unmatched sources (2.6/1.6) as expected, but not 
strongly. The mean off-axis angle of unmatched sources {9 — 6.7 ± 3.9') is indistinguishable from that of matched 
sources (6.2 ± 3.2'; both quoted errors are la RMS). 

For a typical ChaMP source in these fields, the median X-ray (0.5 - 2 keV) fiux is 2.4 x 10~^^ erg cm ~^ s~^, where 
according to the logA^-logS', the cumulative (0.5-2 keV) X-ray source density attains about 30deg~^ fe.g.. lR,osati et af] 
l)2002tl ). The probability of two X-ray sources falling within 10" of each other by chance is < 0.001. On the other 
hand, the number of true multiple sources is of great interest; they may either be lensed images or perhaps binary 

See |\protec1jpEtp: / / asc.harvard.edu/cai /ASPECT/] for a detailed discussion of Chandra aspect solutions. 
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Fig. 3.— Chandra image of the field of ObsID 536 (target MS1137.5+6625). The 16' x 16' AGIST (0.3-5 keV) image (N 
up, E to the left as in all images) was smoothed with a 3 pixel (1.5") gaussian for clarity. 156 X-ray sources detected to 
2 X 10"^'' erg cm~^ s~^ using the ChaMP's X-ray pipeline are marked with circles sized to the Chandra PSF (90% EE). To the 
125 > 2 X-ray sources detected in these 4 chips, we have matched 80 optical sources detected in our combined 4m/Mosaic 
images to r* ~ 25, indicated by 15" boxes. The target cluster is easily visible as the extended source NW of center, and is 
detected as a single source by XPIPE, but like all target objects is excluded from ChaMP analysis. 



Table 4. Source Matching 



ObsID TVx" Nxo"" fxo'' A^xos" fxas"" 
% % 



536 


140 


77 


55 


33 


70 


541 


65 


45 


69 


24 


72 


624 


58 


45 


78 


19 


74 


861 


80 


57 


71 


14 


59 


914 


88 


65 


74 


20 


50 


928 


55 


42 


76 


15 


50 


Total 


486 


331 


71 


125 


63 



''Number of X-ray sources detected in each Chandra image and also having > 2 in counts as described in Sisl 
''Number N or percent / of X-ray sources with high confidence optical counterparts as described in Sl5l 
'^Number A'^ of optical counterparts with spectroscopic classification. 
"^Percent / of r* < 22 optical counterparts with spectroscopic classification. 



AGN, possibly interaction-triggered (e.g., Green et al. 2002). Since these latter types are expected to be rare, a 
large-area survey like the ChaMP is needed to find them. New close pairs of X-ray sources with optical counterparts 
are described in Appendix A. 
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3; Optical, FieldID = .5300 3; x-ray, FieldID = 3300 




Fig. 4. — An example from ObsID 914 (target CL J0542. 8-4100) of the pages used for initial visual inspection of X-ray source 
matches. RIGHT PANELS: A 1' field centered on a detected X-ray source is shown, with N up and E to the left. The 0.3- 
8 keV X-ray image is smoothed with a 2 pixel (1") gaussian for clarity, which facilitates source inspection and examination of 
background and confusion efi'ects. Green circles around the X-ray source have size equal to the PSF (50% EE) at the source 
position. Crosses show the wavdetect X-ray centroid. If other X-ray sources are detected within the field, they are marked by a 
20" box. LEFT PANELS: The corresponding optical r* images shows all optical sources detected by SExtractor, marked with 
3" circles. The Chandra 50% EE circle is reproduced at the X-ray source position. When our automated matching procedure 
finds an optical counterpart, the optical source is marked with a 3" box. X-ray and optical photometric data are displayed to 
the far right. 



6. OPTICAL SPECTROSCOPY 

The nature of many individual sources, especially those with unusual properties, cannot be reliably determined until 
spectroscopy is obtained. Our spectroscopy targets all objects in a representative subset of 40 ChaMP Cyclel-|-2 fields, 
yielding ~ 1000 objects for which we will obtain the most detailed information. With our spectroscopic subsample, we 
will classify sources by type, determine redshifts and hence luminosities and look-back times for evolutionary studies, 
and test photometric determinations of redshift and classification. 

The mean magnitude of optical counterparts is r' ~ 22 for a 23 ksec high galactic latitude Chandra exposure, and so 
fields of this Chandra exposure or larger typically have at least 15 sources suitable for fiber spectroscopy. Multi- fiber 
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Table 5. Spectroscopic Facilities 



Telescope 


Instrument 


Mode 


Grating/ Grism 


A range 


R 


Spectral 










(A) 


(A/AA) 


Resolution (A) 


WIYN^ 


HYDRA/RBS 


multi-fiber 


316@7.0 


4500-9000 


950 


7.8 


CTIO Blanco 4m 


HYDRA 


multi-fiber 


KPGL3 


4600-7400 


1300 


4.6 


MMT 


Blue Channnel 


single slit 


300 1/mm 


3500-8300 


800 


8.8 


Magellan 


Boiler Sz Chivens 


single slit 


1200 1/mm 


4200-5800 


2083 


2.4 


Magellan 


LDSS-2 


single slit 


med/red 


4000-8500^ 


520 


13.5 


Keck3 


LRIS" 


multi slit 


300 1/mm 


4000-9000^ 






FLWO 1.5m 


FAST 


single slit 


300 1/mm 


3600-7500 


850 


5.9 



The WIYN Observatory is a joint facility of the University of Wisconsin Madison, Indiana University, Yale University, and the National 
Optical Astronomy Observatory. 
^Spectral coverage can vary as a function of slit position in the mask. 

•^The W. M. Keck Observatory is operated as a scientific partnership among the California Institute of Technology, the University of 
California, and the National Aeronautics and Space Administration, and was made possible by the generous financial support of the W. 
M. Keck Foundation. 

*LRIS; Oke et al. 1995 



spectroscopy, even on 4 to 6m telescopes, allows us to consistently obtain source classification and redshifts to r' ^ 20, 
and up to l.Smag fainter for objects with strong emission lines. Our overall spectroscopic observing strategy has been 
to use single-slit (FLW01.5m w/ FAST) spectroscopy for isolated r' < 17 sources in shallow Chandra fields, and in 
our deeper fields, which might otherwise introduce scattered light problems for multi-object spectroscopy. Multi-fiber 
wide-field spectroscopy with the HYDRA spectrographs (;Barden:1994) on the WIYN 3.5m and CTIO 4m telescopes 
enables us to acquire spectra for sources with 17 < r' <21. A subset of faint ChaMP sources (21 < r' <23) in fields 
designated for spectroscopic follow-up are being observed with 6m class telescopes such as Magellan and the MMT. 
Table El lists the technical specifications and configurations for each spectroscopic facility used for results presented 
herein. 

Of the 6 fields studied here, each has been observed with a multi-fiber spectrograph on either the WIYN or CTIO 
4m through 2" fibers (Table E}. An average of 15 spectra per field are acquired with this configuration. We obtained 
multi-slit observations for 17 sources with Keck-I/LRIS or Magellan/LDSS2. Single slit observations of 19 sources 
were acquired with the FLWO 1.5m when optically bright (r* < 17) and MMT or Magellan when faint. We perform 
standard optical spectroscopic reductions within the IRAF environment. For multi-fiber reductions, we use the IRAF 
task dohydra. To optimize sky subtraction, we cross-correlated all dispersion-corrected (wavelength calibrated) sky 
spectra from random sky locations in the field against the object-|-sky spectrum from each object fiber using the task 
xcsao within the external IRAF package rvsao. We then median-combine the 9 sky fibers whose spectra correlate 
most closely with the total object -l-sky spectrum. We found this method to be superior to removal of a simple field- 
averaged sky spectrum for faint objects. Slit spectra were reduced using the IRAF task apall with a background 
signal measured locally for each object within the same slit. The multi-fiber spectra are only nominally flux-calibrated, 
since we use just 1 or 2 standards per night down a single central fiber. One or two standards are used for the single 
slit spectra. 

6.1. Object types & Redshifts 

We classify objects spectroscopically simply as broad line AGN (BLAGN; FWIIM> 1000 km s^^) narrow emission 
line galaxy (NELG; FWHM< 1000 km s^^; some of which may be AGN), absorption line galaxy (ALG), or Star. 
Objects with high (> 20) and no identifying features are denoted as BL Lac candidates, as are those whose 
identified Call 4000A region shows a break contrast of < 0.25 (see § 17.4. 3|l . Our classification confidence level is 
denoted as well (0 for not classifiable, 1 insecure, 2 secure). These are all strictly optical classifications. 

Redshifts are measured for extragalactic objects using the radial velocity package rvsao under the IRAF environment. 
This method dep ends on the optical classification. Quasar spectra are quite similar over a broad range of redshifts 
and luminosities l)Forster et alJl2?)0lh . so that - absent the effects of broad absorption lines (BALs) or a strong Lya 
forest absorption - a single high ^ composite serves well as a redshi ft template. For BLAGN, we use the quasar 
template created from over 2200 quasar spectra collected by the SDSS l|Berk et alJl200l)) . We run xcsao over the full 
redshift range in 0.1 redshift intervals, choosing the strongest correlation as the true redshift. Since the correlation 
strength depends strongly on and artifacts such as poorly subtracted sky lines, we inspect and visually verify 
all redshifts and classifications by comparison with overplotted templates. For ALG we use xcsao with a synthetic 
absorption line template. Whenever narrow emission lines are present, we elect to measure the redshift from these 
features to reduce the associated error. For NELG, the task emsao is used to identify emission lines using a line list. 
For stars, we assume an effective radial velocity of zero and allow xcsao to correlate against stellar template spectra 
from ljacobv . Hunter, & Christian ( 1984). Our comparisons between by-eye classifications and these automated results 
suggest an uncertainty of about one full spectral type. 
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Table 6. Optical Spectroscopy 



ObsID 


Telescope 


Instrument 


UT Date 


# of spectra 


536 


Keck 


LRIS 


1 ^ Mav 2nnn 


5 




WIYN 


HYDR A 


n? Anr 2nni 


7 




WIYN 


FTYDR A 


•^0 Tan 900"? 


13 




WIYN 


FTYDR A 


O L U ctll 


4 


541 


WIYN 


HYDR A 


09 Anr 9nni 


12 




WIYN 


FTYDR A 

LL 1 LJ L\t\. 


07 Anr 9001 


9 




MMT 


i_} 1 Lit; v_y 11 CLll 111.^1 


ofi Mav 2nni 

^ U LvL CLy ^\J\J L 


I 




MMT 


BluG cliSyiinGl 


12 July 2002 


2 


624 




HYDR A 

LL L Lj L V Ji- 


1 fi Oct 2001 

LVJ v_yL-L ^Kjyj L 


11 




Magellan 


LDSS-2 


03 Dec 2002 


5 


861 


WIYN 


HYDRA 


19 Sep 2001 


10 




WIYN 


HYDRA 


20 Sep 2001 


2 




MMT 


Blue channel 


11 July 2002 


1 




FLWO 1.5m 


FAST 


15 Oct 2001 


2 


914 


CTIO 4m 


HYDRA 


17 Oct 2001 


15 




Magellan 


LDSS-2 


01 Dec 2002 


5 


928 


CTIO 4m 


HYDRA 


15 Oct 2001 


7 




Magellan 


LDSS-2 


21 May 2001 


1 




Magellan 


BSC 


15 July 2002 


5 




Magellan 


BSC 


15 July 2002 


1 



For X-ray sources with r* < 22 optical counterparts, we have aheady achieved 54% completeness. All sources with 
r > 21 that we have identified spectroscopically to date are BLAGN or NELG, because their strong emission lines 
facilitate classification and redshift determination. We are currently extending to 90% completeness at r* = 21 in at 
least 20 ChaMP fields specifically to map the space density of X-S AGN. 

7. RESULTS 
7.1. Fluxes, Redshifts, and Object Types 

Figure [S] shows a plot of the r* band magnitude against log fx calculated in the 0.5-2 keV band. Overplotted are 
3 lines of constant ^ at 10, 1, and 0.1. We calculate this X-ray to optical flux ratio similarly to lHornschemeier et alJ 
l|2001(l or lManners et all l)2002() . using our SDSS r* magnitudes: 

log^ =log/x + 5.67 + 0.4 r* 

Jr 

and use this as our measure of -P^ throughout. For comparison, we also plot a line showing the typical relation 

Jopt 

expected for 0-S quasars (|Green et al.lll995l) . assuming a constant Uox = 1-5, which corresponds to log = —0.57. 
ChaMP quasars clearly differ from 0-S quasars in the expected sense - X-S objects have typically stronger X-ray 
emission relative to optical (log — —0.15). Alternatively, even optically 'dull' quasars are still detected by the 
ChaMP, since luminous X-ray emission is the primary signature of AGN activity. BLAGN are X-ray bright, being 
largely limited to > 0.1. Galaxies tend to appear towards fainter X-ray fluxes but span a wide range of optical 

brightness. Stars are essentially limited to < 0.1 and r* < 20, which corresponds to a typical M dwarf at distances 
of 400- 1000 pc. 

Apparent trends between fluxes or colors in any survey can be due to a combination of survey sensitivity limits. To 
investigate selection effects in a first approximation, we have generated a very deep simulated SDSS quasar catalog as 
described in LFan. X.lll999j) . based on their observed optical luminosity function and its evolution, and including the 
effects of emission lines and intervening absorption. We have extended the simulation to include AGN with Mb < — 21 
out to a redshift of 6, corresponding to z' magnitudes to 28. We then assume that the mean observed ratio log ^ has 

a constant value of -0.14 with a dispersion a^r = 0.4. (These are the actual values we find for all 41 quasars in our 
sample in the flux regime /x > 5 x 10"^*^ erg cm~^ s~^ where our optical matching rate is high.) From the r' mag for 
each simulated quasar, we then generate X-ray fluxes, luminosities, etc. Finally, we can cut the sample at any desired 
mag or X-ray flux. Primarily, we use a cut of r' < 25 and log fx — —15, similar to our actual flnal sample limits. 
For most figures which require redshifts, we use a cut of r' < 22 more appropriate to the spectroscopic sample. The 
simulated sample is thus based on a quasar OLF with an assumed constant observed y^ for an X-S sample, rather than 
on the actual quasar XLF. Indeed, a prime goal of the ChaMP is to measure the XLF of X-S quasars with improved 
statistics at fainter fluxes where it is not yet well-characterized, and to compare the XLF and OLF particularly at 
redshifts beyond two. 

In Figure |S1 we plot the mean r' vs. log fx values of quasars in this simulation subsample, derived in redshift bins 
of Az = 0.5. The track matches well with the overall trend of BLAGN. There is large dispersion in both optical and 
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Fig. 5. — SDSS r* band magnitudes vs. Chandra log fx (calculated for the 0.5-2 keV band from the broadband 0.3 - 8keV 
counts) for all X-ray sources in the 6 fields studied here. Overplotted are 3 dashed lines of constant at 0.1, 1, and 10 
calculated as described in § 17.11 A fourth (dotted) line marked OSQSOs shows the typical mean values for optically-selected 
quasars (log — —0.57). Large symbols for different optical source types are shown as determined from our spectroscopy: 
broad line AGN (BLAGN), narrow emission line galaxies (NELG), absorption line galaxies (ALG) and stars. Any X-ray source 
not detected in our r*-band images is marked with an arrow at the position of its r* magnitude limit. Small filled circles indicate 
sources with an optical counterpart in at least one filter. The curve traces the mean expected quasar track for z = 0.5 redshift 
bins for the simulation described in §0 imposed flux limits of log /x = — 15 and r* = 25. The bright end point has J = 0.37, 
while the faint end has z — 4.6. Many of the optically faint objects are likely to be high— z AGN. Magnitudes of point sources 
brighter than 18 suffer from saturation effects. 



X-ray fluxes about this mean track. Substantial incompleteness in optical IDs sets in fainter than r* ~ 23.5, where a 
decreasing fraction of X-ray sources have optical counterparts. However, as discussed below in S I7.2I the characteristics 
of X-S quasars clearly differ from the 0-S quasars on which the OLF is based. 

The proportional representation of different object types changes strongly with redshift. Figure shows the overall 
redshift distribution of objects classified in our spectroscopy. BLAGN are from z ^ 0.2 to beyond z ~ 3. We can 
detect and classify quasars to r* magnitudes considerably fainter than galaxies due to their strong broad emission 
lines. We only see galaxies (both NELGs and ALGs) to z ~ 0.8. This effect is due primarily to magnitude limits; our 
spectroscopic limit so far is r* ^ 21 for objects with no strong emission lines within the covered optical wavelength 
range. Figure [3 highlights this by showing that the faintest classified ALGs are at r* = 21. 
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Fig. 6. — Redshift histogram of classified extragalactic sources in these 6 fields. Most objects showing no broad line component 
are below z = 0.8, likely caused by our typical spectroscopic magnitude limit near 21. 



From galaxies in the Hubble Deep Field with good completeness to i? = 23 ijCohen et al.ll200'oj) . the median redshift 
of galaxies at i? = 21 is about 0.5. By contrast, our median r* mag for z = 0.5 galaxies is significantly brighter at that 
redshift. This may not be due to incompleteness of optical foUowup; given both the possibility of an AGN contribution 
and the known correlation between bulge luminosity and black hole mass, we note that X-S galaxies are likely to be 
atypically luminous in the optical band. 

7.2. Luminosities and X-ray /Optical Ratio 

Figure IHl shows the 0.5-2 keV X-ray luminosity of spectroscopically classified galaxies and AGN as a function of 
redshift. BLAGN span a higher range of luminosities than do galaxies, but with som e overlap. Overpl otted is the 
sensitivity of the ROSAT Deep Survey assuming their quoted limit = 5.5 x 10~^^; ijLehmann et al. 1 12001). The 
ChaMP samples luminosities 5-10 times fainter at a given redshift even assuming a more optimistic limit (from their 
Figure 3) of 2 x lO'^^. 

The slope of a hypothetical power law connecting 2500A and 2 keV is defined as aox — 0.384 log-^^ (where these 

are monochromatic luminosities at 2500A and 2 keV), so that (Xox is larger for objects with stronger optical emission 
relative to X-ray. We note that our definition of aox here includes no mean fc-corrections for emission lines, intervening 
or intrinsic absorbers, so aox is directly related to log , which we prefer. 
0-S AGN samples show typical aox values of about 1.5 ijGreen et al.lll995j) corresponding to log ^ = —0.57. We 

For convenience, the exact conversion between the two is cxox = -0.3791og if- + 1.286. 
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Fig. 7. — SDSS r* band magnitudes vs. redshift. Galaxies of all types are seen out to z ~ 0.8, beyond which they are typically 
too faint for spectroscopic classification. BLAGN can be detected at all redshifts and extend well beyond our nominal r* = 21 
spectroscopic limit due to their strong broad emission lines. The wide range of r* mags at each redshift reflects both the breadth 
of object types detected, and the range of exposure time across and between fields. 



find for the current sample of ChaMP BLAGN a somewhat lower value of aox = 1-36 ± 0.15 (log = —0.19). 
iLa Franca et al.l l)1995j) found a similar mean aox of 1.3 for X-S AGN (log ^ = 0). X -ray surveys select X-ray bright 
objects by definition, and so are complementary to 0-S samples. 

A significant correlation between optical luminosity and aox has been noted in many studies based on 0-S quasars 
(Vignali et al. 2001, 2003; Yuan et al. 1998; Green et al. 1995; Wilkes et al. 1994). One cause might be low-energy 
X-ray absorption as reported in quasars at higher redshifts (Vignali et al. 2001, Yuan et al. 2000; Fiore et al. 1998; 
Elvis et al. 1994). Hypothesized to be due to the presence of large gas quantities in primeval galaxies, such absorption 
might also be related to strong dense gas outflows akin to those seen in broad absorption line (BAL) QSOs. However, 
the trend in these samples is primarily with Lopt and not redshift. Another cause of the apparent trend may be a 
combination of wider intrinsic population dispersions in Lopt than in Lx along with superposed selection effects (Yuan 
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Fig. 8. — Wide coverage of the log Lx vs. redshift plane highlights the advantages of varying depths and wide sky area of 
the ChaMP. BLAGN span a higher range of Lx than galaxies but with some overlap. The dashed line shows the minimum 
luminosity corresponding to which the ROSAT Deep Survey was sensitive (see 17.211 . 



et al. 1998, La Franca et al. 1995). 

We detect no discernible dependence of aox on either redshift l'Mathuil l2002D or optical luminosity. Figure shows 
no correlation between log and optical luminosity (which we calculate as vLi, at 2500A, assuming optical spectral 

slope a = —0.5). Galaxies span a lower but equally wide range of log ^ as do BLAGN. The curve in Figure El traces 
the simulated mean quasar track with Az — 0.5 redshift bins and imposed limits of log fx = —15 and r* < 22. All 
input simulated quasars have constant log = —0.12 {aox ~ 1-4). The high luminosity point has z = 0.37, while the 
faint end has z = 4.2. No trend is seen in BLAGN for either data or simulation, because the accretion dominates the 
spectral energy distribution. By contrast, Figure [HI shows a strong correlation between log and X-ray luminosity 
for galaxies. The trend is likely a combintation of 2 effects, one intrinsic, and one extrinsic: (1) The wide range 
of log for galaxies may reflect the unlinking of nuclear X-ray luminosity from the optical luminosity, where the 
latter is strongly affected by host galaxy contributions. (2) The current optical spectroscopic magnitude limit excludes 
sources that are very weak in optical relative to X-ray emission (have large y^) and have small fx- 

7.3. Colors 

The optical colors of ChaMP sources depend on object type, redshift, and reddening. The SDSS provides excellent 
statistics on the colors of the objects that fall into their spectroscopic sample. The ability to compare the photometric 
and color properties of X-S sources with SDSS sources is an important ChaMP feature which will be treated in 
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Fig. 9. — LEFT: log vs. log Lopt shows that galaxies span a range of X-ray to optical flux ratio at least as wide as do 
BLAGN. The curve traces the mean expected quasar track for z — 0.5 redshift bins for the simulation described in § Q and 
imposed limits of log fx = —15 and r* < 22. All input simulated quasars have constant log = —0.3 (corresponding to 
cXox ~ 1.4). The high luminosity point has z = 0.37, while the faint end has J = 4.2. RIGHT: log vs. log Lx correlates 
strongly for galaxies but not BLAGN. The dashed line corresponding to r* = 22 and z — 0.5 is shown for a range of log fx from 
-16 to -12 to illustrate the effect of our spectroscopic limit. The wide range of log for galaxies may reflect the unlinking of 
nuclear X-ray luminosity from the optical luminosity, because the latter is strongly affected by host galaxy contributions. 



upcoming papers. 

Figure E| shows the optical color distribution of ChaMP sources. The mean stellar locus in these SDSS colors is 
shown as two light short-dashed lines. The range of BLAGN is similar to 0-S quasars, where the expected span of z < 3 
quasars in both {g* — r*) and (r* —i*) is from about -0.5 to 1 ((Richards ct al. 2001). We plot the curve of the simulated 
mean quasar track for z — 0.5 redshift bins for our simulation (with imposed flux limits of log fx = —15 and r* — 25). 
The first point on the track to leave the stellar locus is at z = 3.2. Without X-ray information, normal quasars are 
virtually indistinguishable from stars in this color plane until higher redshifts, where quasars redden significantly in 
{g* — r*), with less change in (r* — i*). The reddest {g* — r*) point on the track corresponds to J = 4.6. Galaxies 
are generally redder than BLAGN. For ALGs, the 4000A break is the stronges t color determinant. The mean trac k 
of an E/SO galaxy is shown as a red long-dashed line, using the HyperZ code of iBolzonella. Miralles. fc Pelldl l)2000() . 
assuming a Bruzual-Charlot E/SO galaxy spectrum and typical E(B-V)= 0.05 for Galactic dereddening. The {g* — r*) 
colors of E galaxies redden to {g* — r*) ~ 2 as the break moves redward from the g' to the r' band at z ~ 0.4, and 
the track achieves its reddest {g* — r*) at about z — 1 before looping back blueward to end a,t z — 1.2. Several ALGs 
with 0.4 < z < 0.6 are shown with (g* — r*) close to 3; with log Lx ~ 42, these are likely to be dust-reddened, low 
luminosity AGN. 

The distribution of optical colors changes with X-ray flux as shown in Figure 1111 In our spectroscopic sample to 
date, BLAGN inhabit the full range of X-ray fluxes, but a narrow range of {g* — r*). The brightest X-ray sources are 
all blue in {g* —r*). The curve in FigurelTTl traces the simulated quasar track with imposed flux limits of log /x = — 15 
and r* = 25. The X-ray bright end point has J = 0.37, while the faint end has z = 4.6 Many of the faint red objects are 
likely to be high redshift and/or highly absorbed AGN. X-ray faint blue objects may be the wings of the low-redshift 

quasar ^ distribution. 

Figure IT^ shows SDSS (g' — r') colors plotted against log ^ in the 0.5-2 keV band. Stars are all at log y^ < — 1 
and {g* — r*) < 1.6. Bright BLAGN cluster strongly at high log y^ and blue {g' — r') colors. Objects in this part 
of the plot are largely consistent with 0-S quasars and can be classified as such with high confidence even without 
spectroscopy; a suggested demarcation is shown as the line {g* — r*) < 0.3 log y^ -I- 0.7 in Figure [T^ Redder objects 

with log y^ > —1 but {g* — r*) > 1 (many undetected in our g*-band imaging here) could be one of 2 types. The 
simulated quasar track shows that some are high redshift (z > 4) quasars, but many will be obscured AGN, as shown 
by the high fraction of objects without broad lines in this region. 



7.4. X-ray Bright Galaxies with no Emission Lines 
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Fig. 10. — SDSS ((?*— r*) vs. (r*—i*) color-color diagram for objects in 6 ChaMP fields. Optical objects with no X-ray detection 
are shown as small black points. The typical mean locus of SDSS objects that we use for a sanity check of our photometric 
cali brations is plotted as 2 short-dashed straight lines. The reddest ]3LAGN in both colors is at 2; = 4.93 (Silverma n et al.l 
|2TO2). The solid curve traces the mean expected quasar track for z = 0.5 redshift bins for the simulation described in §|7| with 
imposed flux limits of log /x = — 15 and r* = 25. The first point on the track that is away from the stellar locus is at 2: = 3.2, 
while the reddest (jf — r*) point on the track corresponds to z = 4.6. Without X-ray information, normal quasars are virtually 
indistinguishable from stars in this color piano u ntil z > 3.5 . The mean track of an E/SO galaxy is shown as a red long-dashed 
line, using the HyperZ code of ..Bolzonella. Miralles. fc Pellcl l|200nf) . The (jf — r*) colors of E galaxies redden to (jf — r*) ~ 2 
as the break moves redwar d from the <?' to the r' band at z ~ 0.4. The loop back blueward in (g* — r*) occurs at 2 ~ 1, where 
an L* E/SO has r* = 23.5 JBlanton et aljl200ll'l . 



A small number of X-ray bright, optical normal galaxies ijComastri et alJf2002D have also been found in XMM and 
Chandra surveys. These objects show at best weak optical emission lines. We propose three possible explanations for 
X-ray luminous objects showing no broad optical emission lines. (1) A 'buried AGN' - a nucleus that is intrinsically 
similar to that of a quasar but whose optical continuum and broad emission line region are shrouded by a large 
obscuring column. (2) A low luminosity AGN (LLAGN) where the host galaxy light dominates hght in the spectral 
aperture. One such type of LLAGN may host a nucleus with a massive black hole accreting at very low rates, i.e., an 
advection dominated accretion flow (ADAF). (3) A BL Lac object. 

7.4.1. Buried AGN 

The unified model explains the Seyfert classes 1 & 2 as different views of the same phenomenon l| Antonuccll 1 993j) . 
Quasars are thought to be luminous versions of Seyfert galaxies, so this unified (orientation + obscuration) scheme 
should hold for quasars, i.e., they should have a similar nucleus and dusty torus system in their center. Seyfert 2 
galaxies outnumber Seyfert 1 galaxies 4:1 in the local universe, which is consistent with the covering factor of the 
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Fig. 11. — SDSS {g' — r') colors derived from our NOAO Mosaic imaging vs. Chandra log fx calculated in the 0.5-2 keV 
band. Colors include limits (denoted by arrows) when the object is only detected in a single optical filter. The brightest X-ray 
sources are blue in {g* — r"). The curve traces the mean expected quasar track for z = 0.5 redshift bins for the simulation 
described in §|7| with imposed flux limits of log /x = — 15 and r* = 25. The bright end point has 1 = 0.37, while the faint end 
has z = 4.6 



obscuring tori ("Maiolin o fc R,iek3ll99!TD . and if that ratio can be extrapolated to higher luminosities, then we should 
expect many Type 2 quasars to be detectable in the high-redshift universe. Only a few Type 2 quasars have been 
convincingly detected to date (Akiyama et al. 2002, Norm an et al. 2002, Stern et al. 2001, Sc hmidt et al. 2002), 
implying a strong, perhaps luminosity-dependent evolution ijFranceschini. Braito. & Faddall200 2l). However, infrared 
surveys are currently finding bona fide Type 2 quasars. Several of the IRAS Hyperluminous Galaxies (Beichman et 
al. 1986; Low et al. 1988, 1989) have been shown to be Type 2 quasars (e.g.. Wills et al. 1992; Hines & Wills 1995; 
Hines et al. 1995). More recently, samples selected from the 2MASS include a large number of Type 2 objects with 
quasar-like n ear-IR luminosities. Fully 1/3 of the spectroscopically confirmed AGN in the 2MASS are Type 2 objects 
l)Cutrill20'oH) . Extrapolating to the entire sky, 2MASS s hould detect rough ly 5,000 such objects with K < 15, many 
of which will have IR luminosities in the quasar range llSmith et alJl2002D . Near-IR surveys are less biased against 
absorption and less sensitive to orientation than optical surveys, but will preferentially select samples with high star 
formation rates and copious dust. 

The X-ray characteristics of 2MASS quasars suggest strong absorption (Nh ~ 10^^ cm~^; Wilkes et al. 2002) with 
contributions from direct and unabsorbed, scattered, and/or extended emission. The possibility remains that they are 
intrinsically X-ray- weak. The sample of 2MASS quasars studied with Chandra by Wilkes ct al. (200^ are all at z < 0.3 
(due to the 2MASS K ~ 14.5 limit). At higher redshifts, the reduced effective column in the X-ray bandpass should 
be effective in discovering Type 2 quasars, and help resolve the question of their intrinsic X-ray luminosities. How 
can and should a Type 2 quasar be recognized in an X-ray survey? Not simply on the basis of X-ray hardness. High 
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Fig. 12. — SDSS {g —r') colors derived for Chandra sources from our Mosaic imaging plotted against log in the 0.5-2 keV 
band. The solid curve traces the mean expected quasar track for z — 0.5 redshift bins for the simulation described in §|7| with 
imposed flux limits of log fx ~ —15 and r* = 25. The bright end point has z = 0.37, while the faint end has z = 4.6. All 
input simulated quasars have constant log = —0.3, so the shift towards larger log arises only from the imposed flux 
limits. Objects in the lower right part of the plot are largely consistent with O-S quasars and can be classified as such with 
high confidence even without spectroscopy; a possible demarcation is shown as the dotted line {g* — r*) < 0.3 log + 0.7. 
Clustered with stars at low log y^ are galaxies that likely have little X-ray contribution from an active nucleus. 



column density X-ray-absorbing gas (e.g., a warm absorber) may be interior to the putative obscuring torus, so that 
even normal optical Type 1 (broad line) quasars may show low luminosity an d/or high colu mn abso rption in the X-rays 
(how ever in these cases absorption in the restframe UV is usually apparent: iGreen et al.i r2001: Br andt. Laor. fc Wiiy 
1200 01). Luminosity criteria may be subtle, since the claim for high-L could be based on either Lopt or Lx, perhaps 
including large absorption-corrections. We find objects in our sample that show signs of X-ray or optical absorption, 
or both, but with no strong evidence that these properties are coupled. 

If an X-S AGN is found showing no broad emission lines, oftentimes the search for broad lines may be very limited, 
including only what is seen on a discovery optical spectrum. A search for broad Ha is somet imes possible from the 
ground, and may reveal instead a luminous Type 1.9 quasar ijAkivama. Ueda. fc Ohta|l2002f) . but the classification 
is by nature somewhat arbitrary. Spectropolarimetry can reveal broad line flux, but since those same photons must 
also be in the total flux spectrum, broad line detection becomes a question of adequate . Indeed, many IR-selected 
Type 2 quasa rs show lower frac tional polarization than Type 1 quasars, presumably because of strong stellar light 
contributions ijSmith et al.ll2003) . Furthermore, there are few objects for which we can convincingly claim an optical 
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Type 2 classification without examining line ratios in detail in spectra of higher ^ . Instead we define a "buried AGN" 
here as an object that has either no or only narrow emission lines in optical spectra, strong evidence for logNn > 22 in 
the rest-frame, and log Lx > 43 without absorption-correction. If the signature of X-ray absorption is clearly detected, 
then an absorption-corrected log Lx > 44 qualifies as a Type 2 quasar candidate, with confirmation pending high jj 
foUowup spectroscopy. We find no such objects in these 6 fields. 

7.4.2. Low- Luminosity AGN 

Many of these X-ray selected galaxies are so distant that their angular diameters are comparable to the slit (or 
fiber) widths used in gro und-based spectroscopic observations. U sing integrated spectra of a sample of nearby, well- 
studied Seyfert 2 galaxies. iMoran. Fi lippenko. fc Chornocld l)2002f) demonstrate that the defining spectral signatures of 
Seyfert 2s can be hidden by light from their host galaxies. At z ~ 0.6, a 1.5" slit encompasses a region about lOkpc, 
which includes a substantial fraction of the host starlight emission. Some 60% of the observed objects would not be 
classified as Seyfert 2s on the basis of their integ rated spectra, which is con i parable to the fract ion of Chandra sources 
identified as "normal galaxies" in deep surveys l|Hornschemeier et alJl200'll iBarger et alJl2002D. 

These 'buried AGN' may host low rate (advection dominated) accretion flows fADAFs: iNaravan fc YilllQQIifl . The 
X-ray spectra produc ed by ADA Fs are relatively hard, so ADAF sources may contribute a fair share of the hard 
(>2keV) background l)Yi fc Bou ghn 1998). Half of the 2-10 keV CXRB could consist of low-luminosity (Lx < lO^'^ 
ergs/s) sources if the comoving number density is ~ 3 x 10~'^Mpc~^, comparable to the density of L* galaxies (e.g., 
Peebles 1993). ADAF sources should be characterized by inverted radio spectra Ii, oc 

There are 8 optical counterparts to X-ray sources in the current sample that have hard band^^ luminosities logi(2- 
8 keV) > 43 and no signs of broad optical emission lines in their spectra. The luminosity in this band is plotted against 
hardness ratio in the right panel of Figurell3l (A similar Lx criterion in the 0.5-2 keV band would be logL(0.5-2 keV) 
> 42.4) Half of these Type 2 candidates are hard, with HRq > 0. Optical spectra of this sample are shown in FigurelT^ 
Several of these objects have spectra of absorption line galaxies with only a weak narrow [OH] emission line. Several 
counterparts are also classified as narrow emission line galaxies with strong lines of [OH], [OIII] and H/3. All these 
objects are at redshifts less than unity. 
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Fig. 13. — Hardness ratio vs. X-ray luminosity in 2 bands. Most of the luminous objects (quasars) are soft. This may arise 
both because they are at high redshift (where absorption causes smaller changes in hardness) and because they are intrinsically 
less obscured. RIGHT: Plotting against a harder LxC^-S keV) does not significantly change the plot. A dashed line shows our 
criterion for objects with no broad lines to be classified as 'buried AGN'; logLjf (2-8 keV)= 43. 



Newer CXRB models attempt to include the observation that most hard faint sources from Chandra and XMM 
Newton surveys are objects at z < 1, which requires a very different evolution rate for Type 1 and Type 2 objects. 
This must at least modify the meaning of Unification; all Type 2 AGN are not necessarily Type 1 from a different 
angle. The covering fraction of the putative obscuring torus may also change, and probably evolves statistically with 
cosmic time. The relative dearth of high redshift Type 2 quasars could be explained by selection effects. As discussed 
in S 17.11 these objects may simply fall below our optical detection limits, or at least our r* = 22 spectroscopic limit. 
Since these objects have low surface densities, as the ChaMP area grows, a volume-limited subsample will help define 
their relative space densities. 

While our H band is nominally 2.5-8 keV, we use 2-8 keV luminosities here to facilitate comparison to other studies. 
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Fig. 14. — Spectra of ChaMP candidate absorbed AGN. These optical spectra correspond to 4 X-ray sources with luminosities (logLjf (2- 
8 keV)> 43 erg s"!), hard X -ray spectra, and no signs of broad emission lines characteristic of optically-selected quasars. Grey regions 
mark the telluric absorption bands. 



7.4.3. BL Lacs? 

The first known BL Lacs provided the three observational characteristics of the class (e.g., 
iStein Odell. fc Strittmatterlll97()t lAngel fc StockniaiJll98(]t) . repeated here as presented bv lJanniizi. Smith fc Elstonl 
l) 199.11) : "1.) They are intrinsically luminous with strong and 'rapid' variability. (BL Lac objects have been observed 
to vary at X-ray, optical, infrared, and radio wavelengths.) 2.) They have 'featureless' optical spectra. (How to define 
"featureless" has been a question for debate [...]). 3.) The electro-magnetic radiation from BL Lac objects is observed 
to be linearly polarized." In general the polarization requirement has been for the observed degree of polarization to 
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be larger than that possible from differential extinction (at most wavelengths a percentage polarization of greater 
than 3%) or for the degree and/or position angle of the polarization to be variable. Together, these observationally 
based criteria reflect what we now know is the distinguishing physical characteristic of all BL Lacs, the presence of a 
strong source of synchrotron radiation in these sources. 

We are interested in the identification of BL Lacs because it is possible that with only single epoch photometry and 
spectroscopy (and no polarimetry data), objects which we have currently identified as buried AGN might be BL Lac 
objects. Lacking strong emission li nes, the optical and nea r-IR spectra of many BL Lacs (particularly X-ray selected 
objects; see EMSS survey results of lMorris et alJlT991l and lJannuzi! Elston fc Si^lT99l can be dominated by host 
galaxy starlight and look very similar to what we call absorption line galaxies. 

While more data on these objects would be advantageous, we have checked whether BL Lacs might be hidden in 
our absorption line galaxy sample. First, we chec ked for k nown radio emission fro m our X-ray sample. Since all 
known BL Lac objects are strong radio sources (e.g. I.TannuzilllQQ fl: Stoc ke et al.lll99f]() . the presence of strong emission 
(from a compact synchrotron core) might indicate the presence of a BL Lac. Unfortunately, no data are available o f 
sufficient se nsitivity for our AL Gs. Our cross-correlation of our X-ray source fist to the FIRST (Whiteir^ 11993) 
and NVSS l|Condon et all [199^ radio surveys with a 5" search radius did yield 2 matches to other sources. One is 
a bright nearby radio-quiet Sy 1 and the other is likely to be a distant radio- loud quasar (see Appendix^). Second, 
we checked the strength of the 4000A break in our spectra of all sources we classed as ALGs or BL Lacs (all those 
with the strongest emission line equivalent width Wx < 5A). As discussed by the EMSS team, a weak break might 
indicate the presence of non-thermal (featureless) emission contributing to the total emission - again a possible sign 
of a BL Lac being present in the source. We measured the 4000A break contrast using the using the definition of 
Dressier fc ShechtmanI l)1987j) 

Br • = 

4000A / + 

where is the average flux in the rest wavelength range 4050-4250A and /" the average flux between 3750-3950A. 
As long as the in this wavelength region exceeds 5 per pixel, we count break strengths < 0.25 as BL Lac candidates 
pending confirmation from any of the 3 criteria listed above. 

We found just one optical spectroscopic BL Lac candidate; CXOMP J05421. 5-410206, with r* = 19.6 at z = 0.637. 
While in the same range of luminosity as known BL Lac objects, it has a relatively low log — —0.6. We note that 
for optical magnitudes V < 20, the optical l uminosity function for X-ray-selected BL Lacs yields a surface density of 
about 160 steradian"^ or about 0.05 deg~^ (iWolt er et al.lll99^ . The optical QSO surface density is about 600 times 
as high at similar magnitudes (jMever et alji200 1'). so few BL Lacs are expected overall in the ChaMP. 

7.5. Signs of Absorption? 

Figure El shows HR values corrected as described in § O Many previous studies have assumed that when total 
source counts exceed 30 or 50, that the HR values are reliable without accounting for errors. In fact, the combination 
of source spectrum and telescope effective area may require a surprisingly large number of counts just to achieve 
detection in all the bands composing an HR calculation. The fraction of limits (mostly upper limits, since Chandra's 
H band sensitivity is significantly lower) decreases below 10% only above about 120 counts (for these observations 
typically at about log fx = —14 in Figure [T5|l . About half the sources have HRq upper limits even for sources with 
30 counts or more. Analysis of only those points detected in both bands would lead to an apparent hardening at 
low fiuxes. Correctly incorporating the upper limits via survival analysis (ASURV Rev 1.1 Isobe & Feigelson 1990; 
LaValley, Isobe & Feigelson 1992), we find no significant evidence for hardening of the spectra toward lower fiuxes. 
Excluding sources with fewer than 30 (0.3-8 keV) counts, the median (mean) of HRo for sources with < 120 counts is 
-0.66 (—0.58 ± 0. 02), while for sources with more counts we find -0.62 (—0.56 ± 0.03). This is not inconsistent with 
the conclusion in Kim et al.l l|2003bD that spectral hardening at faint fiuxes is significant only in the softest band (SI; 
0.3-0.9 keV) and thus most likely due to absorption. While spectral hardening towards fainter fiuxes may be real, 
careful accounting for errors in hardness ratios is important: spectral stacking (e.g., Brandt et al. 2001; Tozzi et al. 
2001) is a powerful alternative for these analyses. 

Assuming that F ~ 1.9 for most BLAGN, the hardness ratios indicate that most BLAGN likely suffer some absorp- 
tion, whereas strong absorption is rare. A much larger fraction of NELG and ALG show evidence for logNn > 22. 

X-S quasar samples may be less biased against absorbers (both intrinsic and line-of-sight) than are 0-S samples. 
This advantage is expected to be especially important at high redshifts for two reasons. First, X-rays from high-z 
objects can be detected through a larger intrinsic absorbing column of gas and dust because the observed-frame X-ray 
bandpass corresponds to higher e nergy, more penetrating X -rays at the source; the observed-frame, effective absorbing 
column is ~ A^h/(1 + z)"^'^ (jWilman fc Fabian 1999). Second, optical color criteria for high-z candidates (e.g., 
blue filter dropouts) tend to select objects with strong Lya forest decrements. 

Do red optical colors reliably predict hard X-ray spectra or vice- versa? Figure [TBI shows two quantities, HRq and 
(g* — r*), which may serve as crude proxies for absorption measurements. While {g* — r*) is affected by dust via 
optical reddening, it is also strongly dependent on object type, host galaxy fraction, and redshift. While the hardness 
ratio HRq is affected by the total X-ray absorbing column, it is also dependent on redshift since the effective column 
decreases as noted above. We find no overall correlation between these two absorption indicators, as is evident from 
Figure El X-ray and optical absorption are effectively decoupled for several physical reasons: 1) The gas-to-dust 
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Fig. 15. — Corrected hardness ratios HR = [H — S)/[H + S) as a function of X-ray flux. HR values here are corrected to 
represent a source on-axis on a front-illuminated (FI) chip, and to include upper or lower limits to HRo (arrows), based on 
2(7 upper limits to H or S counts, respectively. Horizontal dashed lines show the HRq corresponding to a typical unabsorbed 
AGN (bottom; F = 1.9/Nh ~ cm~^), the average hardness ratio of the total CXRB spectrum (middle; F = 1.4/Nh = or 
F = 1.9/logNH > 21.5), and a highly absorbed AGN (top; F = l.g/logNn > 22). Most BLAGN likely suffer some moderate 
absorption, whereas a large fraction of NELG and ALG show evidence for log Nh > 22. The brightest 8 sources between 
— 13.4 < log fx < — 12 (typically soft) are left outside the plot range for clarity. 



ratio varies, depending on densities, abundances, and ionization level ijMaiolino et alJl200H) . 2) Physic ally d istinct 
absorbing regions tWcingar tner fc Mu rray 20d2|). 3) Variations in geometry, reflection, and self-shielding l)Elvisl l20001 . 

We note that for the 35 X-ray hard objects (HRq > 0) suggestive of strong absorption, there is some weak evidence 
for a trend in Figure El Using ASURV (LaVallcv. Isobe & Feigclson 1992), we include the limits in both axes,^** and 
find no significant correlation (probability of no correlation is 9.8%). At high columns, other scatter-inducing factors 
such as redshift, orientation, and/or reflection may be overcome by the effects of intrinsic dusty absorbers. This 
suggests that the correlation should be tested with a larger sample, as will become rapidly available from the ChaMP. 

Hard red objects in the upper right hand corner of the plot (e.g., HRq > and {g* — r*) > 1.5) likely contain a 
high fraction of obscured AGN; none with spectral classifications to date are BLAGN. Soft red objects in the lower 
right (e.g., HRq < and {g* — r*) > 1.5) could be high redshift AGN (where unabsorbed X-rays are shifted into 
the Chandra passband) or soft galaxies above z 0.3. Spectra of the hard blue BLAGN (e.g., HRq > —0.2 and 
(g* ~ r*) < 0.5) show mostly weak narrow emission lines, and several appear to have significant blueshifted CIV 

^'^ We employ a generalization of Kendall's nonparametric rank correlation test developed by Brown, Hollander & Korwar (1974, hereafter 
BHK test). This method assumes no distribution for censored data and allows either type of censoring simultaneously in both axes. 
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Fig. 16. — SDSS {g' — r') colors derived from our 4m Mosaic imaging vs. corrected Chandra hardness ratio HRo- Since a 
range of redshifts and optical object types may pertain, {g* — r*) is at best a crude measure of optical reddening - as is HRq for 
X-ray absorbing column. Nevertheless, above HRo ~ a weak correlation exists between these indicators (P = 10%; see S I7.51 . 
This suggests that at high columns, the effects of intrinsic dusty absorbers may dominate the dispersion induced by factors such 
as redshift, orientation, and/or reflection. 



absorption. These objects are sim ilar to BALQSOs, which are only slightly reddened l)Reichard et alJl2003j) but known 
to be heavily absorbed in X-rays l|Green et alJl200lt) . 

7.6. Clusters of Galaxies 

Clusters are important probes of cosmology, especially at high mass and redshift. In general, a surface density 
of > 0.003 deg~ ^ for clusters more massive than Coma above z > 1 violates ACDM cosmologies at 95% confidence 
( Evra^^^y|2002). Evolution of and correlations between X-ray temperature, luminosity, and gas mass in clusters 
i Vikhlinin et a l. 2000) enables the strongest available constraints on hierarchical models. Chandra's high sensitivity, 
low background and good spatial resolution facilitate the detection of spatially extended emission which at high galactic 
latitudes corresponds to emission from groups or clusters of galaxies. 
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Our prototype X-ray extended source detection algorithm works by performing photometry in annuh around sources 
detected by wavdetect, and compares the best-fit gaussian and beta models to the expected PSF size. All of the 
targeted clusters are detected. Some clusters also have point-like emission from a cD galaxy or quasar, which can 
confuse the automated algorithm. These are generally noticed during our visual inspection, which also provides an 
opportunity to flag any X-ray source surrounded by a concentration of optical galaxies. Using this approach, we have 
found 3 serendipitous clusters in these 6 fields, each of which described further below. The ChaMP makes feasible 
simultaneous optic al and X-ray searches for clusters. Further refinement and testing of cluster selection is underway 
ijPerlev et al.ll200'^ that combines automated and by-eye X-ray source extent determinations, Voronoi Tesselation and 
Percolation analysis on color-filtered optical images, and red sequence identification in optical color-magnitude plots. 

7.6.1. CXOMP J033912.2-352614 

During our standard visual inspection of each X-ray source, we noted extended X-ray emission centered near 
03''39"12.2^ - 35°26' 15.0". Figure [T7I shows the X-ray and optical r' band image of the region. A model to 
fit the cluster profile is unconstrained, but fixing to an average cluster value oi (5 = 0.6 y.oncs & Forman 198^, we de- 
rive a core-radius = 47±12" and 2821 counts (0.7-2 keV).^^ The (0.5-2 keV) flux is therefore 3 x lO'^^erg cm'^ s'^ 
This estimate may be affected (at the ^ 20% level) by emission from the extended halo of the central bright (10th 
mag) galaxy of the Fornax cluster, NGC 1399, which lies 9' to the west outside the Chandra field of view (Paolillo et 
al. 2002). 




Fig. 17.— CXOMPJ033912. 2-352614. LEFT: We note extended X-ray emission centered near a = 03 39 12.2, 5 = -35 26 15 
in this ACIS-S image from ObsID 624. The field is approximately 2' NS, and 1.3' EW. The orientation (in this and all images 
herein) is N up and E to the left. Circles show the positions of detected X-ray sources, with circle size corresponding to the 
Chandra 95% EE at this off-axis angle. RIGHT; Our CTI04m/Mosaic imaging indeed shows a cluster of galaxies just N of a 
V — 12 star. The bright star is not the X-ray source since it is 19.5" South of the X-ray centroid, and stars do not generally 
produce diffuse X-ray emission. 

In our optical imaging, the source corresponds to what appears to be a cluster of galaxies just N of a bright {V = 12 
mag) star. Figure [T7I shows the optical r' band image of the cluster next to the (2-pixel gaussian) smoothed X-ray 
image. We are confident that the X-ray emission docs not originate from the star. Stars do not produce significantly 



The cluster is only detectable to about 3rc, within which there are 1630 counts. 
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extended X-ray emission. Even if they did, the star could not reasonably be as sociated with the di ffuse emission situated 
as it is at 1.7 core radii (20") from the emission centroid. As explored bv iPerlev et alJ l)2004D . the {g* — r*) ^ 2.5 
colors of galaxies within 18" of the X-ray centroid define a red sequence corresponding to cluster elliptical galaxies at 
z~ 0.4 -0.5. 




Fig. 18.— CXOMP J033841.4-353133. Smoothed ACIS-I (LEFT) and optical r' image (RIGHT) of the field in ObsID 624 
containing a bright background quasar (see SrO)l. Each image is oriented N up, E to the left, covering 6.7'N-S and 4.3'E-W. 
The BLAGN a,t z — 0.3597 is marked by a box in the optical image. Circles represent the Chandra 95% EE radius. The large 
(37") EE radius and apparent non-point morphology of the X-ray source reflect the PSF at this large (12.6') off-axis angle. 



7.6.2. CXOMP JO5424O.9-4O5515 

The bright EXOSAT and ROSAT AU-Sky Survey X-ray source EXO 0541.0-4056 was studied in the ASCA Medium 
Sensitivity Survey as lAXG J054242-4054. Our spectrum reveals a BLAGN at 2; = 0.7243, but the X-ray image clearly 
shows extended emission out to > 40" from the quasar (Figure^)). A different stretch and zoom of the same region is 
shown at the top of Figure^ To better examine the cluster emission, we first created an image in the 0.7-2 keV band 
only, and masked out all point sources and the quasar with a circle of 4.5"radius. After smoothing the masked image 
with a 6" Gaussian, we determined the peak of the diffuse emission to be at 05^42'"40.4* - 40°55' 7.0", the NE of the 
quasar, corresponding to a projected distance of 80kpc, assuming the quasar and cluster are at the same redshift. 
Since the cluster emission centroid is offset by nearly Tc from the brightest cluster galaxy (the quasar), this cluster is 
not relaxed, and is most likely still sloshing fro m a merger. The fracti on of such offsets as a function of redshift is 
potentially a powerful probe of cluster evolution ((Gohn. Bagla. fc Whita,,2001'l . 

Extracting a profile around the cluster centroid, we find the cluster is well-fit with a /9-model of core radius Vc — 
12 ± 2.5", with j3 = 0.6 ± 0.07 and total flux of 340 ± 40 counts^^, fitting down to the background intensity of 
0.03 counts/pixel. The corresponding cluster flux of 3.6 x 10"^'' erg cm~^ s~"'^corresponds to a cluster luminosity of 
lO'^^'erg s"i. 

Deep imaging at fine spatia l resolution also reveal s evidence in the opti cal and near-infrared for clusters of galaxies 
around (radio loud) quasars l)Yee fc ElHngsonlll993D at similar redshifts l|Sanchez fc Gonzalez-Serr;iiioll2002D . These 
quasars often do not inhabit the cluster cores, which may support the interaction/merger hypotheses for their fueling, 
since lower velocities at large radii enhance merging efficiency (Bckki &: Noguchi 1994) . 



The errorbar includes the uncertainty due to the value of the (3. 
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Fig. 19.— CXOMP J054240.9-405515. This source has weU over 5000 counts in this ACIS-I image (LEFT), and shows a 
sharp peak with significant (~ 40") extent (see also the top panel of Figure |1] for a close-up). Our optical spectrum of the 
corres pondi ng r* = 18.3 object (RIGHT) centered in each of these images (230" N-S, 140" E-W) reveals & z = 0.724 BLAGN 
(see ? 17. 6. 21 . The extended emission quite likely arises in a cluster of galaxies whose members are several mag fainter than the 
quasar. 



7.6.3. CXOMP J234815. 7+005350 

The X-ray emission from this object is consistent with a point source, but the optical image looks convincingly like an 
intermediate redshift group or cluster of galaxies (Figure EHll . We tentatively assign a redshift of z = 0.41 to the cluster, 
based on a spectrum (Magellan 6.5m with LDSS2 03 Dec 2002) of the brightest r* = 19.7 galaxy. The source flux of 
fx = 1-0 X 10~^^erg cm~^ s~^ at this redshift corresponds to log Lx = 41.84. If there is no significant contribution 
from cluster emission, the evidence suggests that an X-ray bright elliptical galaxy dominates the group. The object is 
not likely to be a BL Lac because the Call break strength is normal for an elliptical galaxy (0.57 ±0.1). Some diffuse 
optical emission may be from the bright extended disk of the host galaxy (~ 15kpc in radius), in which case a knot of 
optical emission near the X-ray centroid suggests the possibility of an ultraluminous X-ray source (ULX) in the host. 
Howe ver, this would be 30 times brighter than the brightest known ULX in the Antenna galaxies l)Zezas fc Fabbianol 
|2002|) or M82 (Kaarct ct al. 2001J. Also, we would expect to see strong optical emission lines from the associate star 
formation regions in the disk. 

8. MULTIPLE SOURCES 

Close X-ray sources (< 10") with optical counterparts are good lens candidates, especially if there is also extended 
emission signifying a mass ive potential. They are also interesting for study of the effects of interactions/mergers on 
activity l)Green et alJ200^ . Sources with larger separations wh ere one lac ks an optical counterpart are good candidates 
for X-ray jets, whose emission barely dims at all with distance l)Schwartai2002.1 . We have so far identified 2 close pairs 
of X-ray sources in these 6 fields. 

We find a new close pair of X-ray sources with optical counterparts in the field of ObsID 861, as shown in Figure!^ 
A Magellan spectrum from 14 July 2002 reveals that the brighter object CXOMP J234835. 3-1-005832 is & z = 0.95 
BLAGN. Since the 90% encircled energy (EE) PSFs of these sources overlap, we estimated their individual counts as 
follows. We first measure the counts for each source in a core radius of 3pixels. We then measure the total counts in 
the union of 2 aperture of radius of 10.4 pixels (5") encompassing both sources, and apportion the total counts to each 
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Fig. 20.— CXOMP J234815. 7+005350. LEFT: The X-ray source has just 31 ± 8 counts in this ACIS-S image and morphology 
consistent with a point source. Circles show the 95% EE PSF size, which has radius 6" at this off-axis angle (4.8'). RIGHT: 
The r' band image shows a group of galaxies coincident near the X-ray centroid. From our optical spectrum (Magellan 6.5m 
with LDSS2 03 Dec 2002), we tentatively identify the brightest r* — 19.7 galaxy just NW of the X-ray centroid as a z = 0.4 
absorption line galaxy (see S 17. 6. 31 . Both images are 240" N-S, 90" E-W in size. 

according their core counts ratio. The fainter source CXOMP J234835. 5-1-015836 is beyond our spectroscopic survey 
limit, but should be followed up as a lens candidate. The 2 objects' values differ significantly, (with log = 0.72 
and 0.34 for the fainter and brighter object, respectively) so we believe that they are unlikely to be lensed. However, 
not all lensed objects are achromatic across a wide range in frequency; microlensing, time variability and differential 
absorption are possible reasons (e.g.. Green et al. 2002; Morgan et al. 2001). 

CXOMP J234758.0+010329 is an r* = 20.1 ALG at z = 0.2481. As shown in FigureE21 CXOMP J234758.9-)-010343 
and CXOMP J234758. 0-1-010329 are separated by only 10", and both appear to belong to a group of galaxies. We 
detected no extended emission from the group. 

9. PLANS AND PROSPECTS 

Significant (> 10%) cosmic variance occurs between the logA^logS' in imaged X-ray fields ijRosati et alJl200^ at 
brighter fluxes (log fx > —14.3 cgs; 0.5-2 keV), where cumulative source densities are just 130deg~^ (~10 sources per 
ACIS-I field). How strong is the X-ray cosmic variance, and is it associated with cluster- mass lenses? Using Chandra 
data, an X-ray point-source excess over the ex pectations of blank fiel ds in the fields of se veral clusters has been 
claimed f Cavm et al. Il2001t iMartini et al ."2002^. In the field of A2104. iMartini et all l)2002() found an unexpectedly 
high fraction of these sources to be red cluster member galaxies both with and without optical spectroscopic signatures 
of AGN. T he ChaMP has a lar ge fraction^^ of target clusters in the field list, and so is well-suited for testing for any 
such effect. iKim et al] ()2003by) examine ChaMP cluster and non-cluster field subsamples and compare the co-added 
logA^-log^, finding no significant difference in point source counts. 

In these 6 fields we have discovered 3 new clusters of galaxies. At least two are above z = 0.4, while the third 
still lacks a redshift. Two show X-ray extent, of which one hosts a bright quasar. Further spectroscopy is needed 
to verify membership of faint galaxies seen optically. Our optical imaging program will allow a new generation of 



Of the full field list, 25% have cluster targets. Of the fields with 20 ksec exposures or greater, half have cluster targets. 
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Fig. 21.— The pair of X-ray sources shown here - CXMOP J234835.3-I-005832 and CXOMP J234835. 5+005836 - are separated 
by 4.9" in the ACIS-S image (LEFT), and both have rehable optical counterparts in our Mosaic r'band image (RIGHT). The 
field shown is approximately 30" NS, and 18" EW. Our spectrum of the brighter more southern object shows it to be a BLAGN 
aX z = 0.9469. Spectroscopic follow-up of the fainter object i* = 23.3 should help determine whether the system is lensed. 



studies of the clustering of X-S AGN and of AGN-galaxy correlations. With their low surface density and high 
luminosities, do AGN trace high peaks of the density field, or do they randomly sample the galaxy distribution? 
Is there a change in form of the AG N correlation function at small separations, implying interactions, or triggering 
(^fstathiou & Rces 1988; Kocliaiick, Falco fc MunodlToOgj) ? What is the lifetime of activity in galactic nuclei? With 
the assumption that luminous quasars reside in the most massive host halos, weaker clustering implies shorter quasar 
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Fig. 22.— X-ray (LEFT) and optical r'-band (RIGHT) images of CXOMP J234758. 9+010343 and CXOMP J234758.0+010329. 
Circles show the PSF size (50% EE) at this (7') off-axis angle. Each image is 90"NS and 70" EW. The apparent X-ray extent 
of the individual sources is consistent with the Chandra PSF. The southern object CXOMP J234758.0-I-010329 is a r* = 20.1 
ALG at z — 0.2481. The optical counterpart centroids are separated by 19". While these X-ray emitting objects may belong to 
a group or cluster, we detect no significant extended emission. 

lifetimes fflai man.. fc Hui 1120011: iMartini fc Weinberg|l200H) . 

Multi-wavelength data, beyond the optical and X-ray, is also being planned in the longer term. A subset of ChaMP 
sources are too faint optically to be identified in our imaging program, but we have initiated ground-based near-IR 
imaging. This reveal many sources which are optically faint due to obscuration likely including a significant subset 
of the AGN population. We have also begun a program to obtain high resolution radio foUowup observations of 
the bright radio sources (NVSS and/or FIRST identifications) in the ChaMP sample using the VLA A array at 3.5 
cm. This not only allows us to study the radio properties of the X-S sample, but also assists in source classification. 
Most of the sources observed to date are unresolved in the radio, as expected for typical broad-lined AGN. However 
a significant subset have complex double/triple structure and another subset prove to be mis- identifications due to 
the low resolution of the radio survey data. Once these observations are complete, we anticipate both extending the 
wavelength coverage for these bright radio sources to 20 cm and planning deeper observations to study the radio-quieter 
part of the X-ray sample. 

Beyond the mechanics of source identification and classification, a primary science goal of the ChaMP is to obtain an 
AGN sample which is more representative of the AGN population as a whole than previous samples, due to selection in 
a broad, hard X-ray band. With this sample we seek to understand the intrinsic population, study the range of spectral 
energy distributions and fractions of various source types in order to better constrain models for the geometry and 
physics of the inner regions. To some extent such studies will be facilitated by the multi-wavelength data described 
above, but in the longer term we also expect to follow-up one/more well-defined subsets in more detail to better 
understand the emission mechanisms involved in the various wavebands. 

9.1. Public Data 

Quantities discussed and plotted here are available in catalog form along with image data for the 6 fields in this 
paper through our web site \protect ittp://hea- www. harvard.edu/CHAMP/ For scientists using ChaMP catalog or 
image data, we appreciate an acknowledgment: "This work made use of images and/or data products provided by the 
Chandra Multiwavelength Project (ChaMP; Green et al. 2002), supported by NASA. Optical data for the ChaMP 
are obtained in part through the National Optical Astronomy Observatory (NOAO), operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science 
Foundation." 
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(CXC). 

Wc arc indebted to the staffs at Kitt Peak, CTIO, Las Campanas, Keck, and MMT for assistance with optical 
observations. Thanks to Arne Henden for accurate coordinates for standard stars, and to Warren Brown, Perry 
Berhnd, Michael Calkins, and Susan Tokarz for providing or helping with spectroscopy on Mt Hopkins. Thanks to the 
referee for a thoughtful and thorough review. 

Hearty thanks to the whole Chandra team for an exquisite observatory and a quality archive. We are particularly 
grateful to the late Leon VanSpcybrocck for allowing us early access to his proprietary Chandra data, and for the 
astonishing design quality of the Chandra, mirrors. Leon uniquely combined brilliance, integrity and humility. That 
inspiration and his legacy, Chandra, live on. 
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APPENDIX 

INDIVIDUAL OBJECTS IN 6 CHAMP FIELDS 

Below we note idiosyncracies of particular fields and objects therein, if they were not discussed previously. We also 
mention previously known sources, based on a NED search within 12" of all X-ray sources. We do not tabulate or 
reference sources that have been detected in the growing number of all-sky or serendipitous surveys (e.g., ROSAT, 
SDSS, 2MASS, APM) unless a spectroscopic classification has been published. 

MSim. 5+6625 

At z = 0.7824, the tar get cluster MS1137.5+6625 (Leon van Speybroek, PI) is the second most distant of the EMSS 
clusters l|Donahuell999j) . and is part of a GTO program to combine X-ray and radio measurements to cons train Ho and 
Op thr ough the Sunyaev Zel'dovich effect. The Chandra image of the target cluster has been analyzed bv lCrego et alJ 

CXOMP J114148. 0-1-660603, with 125 H and 19 S band counts, this is the hardest amongst all the sources in this 
study {HR = 0.74 ± 0.08). It is relatively bright (r* = 21.9, i* = 20.5), but not detected {g* > 25.3) in the g' band. 

CXOMP J114013.9-f 655940: The very re d color ((q' - r') > 2.9) of this r* = 21 absorption line galaxy is unexpected 
for a 2; = 0.5 elliptical galaxy such as this l)EisensteirJl2001j) . and even more so given that the object has a very soft 
X-ray spectrum {HRo < —0.67). With logLx = 42.56, this galaxy certainly contains an AGN, but the optical light 
may be dominated by a reddened host galaxy. 

V14 16+4446 

V1416-I-4446 is a z = 0.40 cluster included in the same GTO program (Leon van Speybroek, PI) as MS1137-f 6625. 

Included within the ACIS-I field is the bright z = 0.114 Seyfert 1 galaxy PG1415-f451. Here there are 5067±72 
ACIS counts, corresponding to 9 x 10~^^erg cm~ ^ s~^ (0.5-2keV). T his source is also detected in the Faint Images 
of the Radio Sky at Twenty cm (FIRST) survey ijWhite et alJll99''^ . with 1 ± 0.14 mJy, but is nevertheless radio- 
quiet (|Xu. Livio. fc Baumlll999|) . Using a 5" search radius, the only other FIRST source we find in these 6 fields is 
CXOMP J141655.7-h445452. With r' = 20.94, an (0.5-2keV) flux 2 x lO-^-^erg cm'^ s'^ and 1.37 ± 0.34 mJy in the 
FIRST survey, it is likely to be a radio-loud quasar. 

LP944-20 

While the target of this observation (PI Gibor Basri; sec Rutledge et al. 2000) is a brown dwarf star, the field is near 
(~ 13' from) the center of the Fornax cluster, which includes the two bright galaxies NGC 1399 and NGC 1404. The 
field has been studied in the optical in detail recently by Hilker et al. (1999a,b) who discovered a background galaxy 
cluster at z=0.11 behind the center of the Fornax cluster as part of the 2dF Fornax Cluster Spectroscopic Survey 
(FCSS; Drinkwater et al. 2000). 

FCSS J033922. 5-352530 (CXOMP J033922.4-352530) is a = 17, r* = 17.56 galaxy at z = 0.04664. 

NGC 1404 a,t V = 1947 km is an 11th mag elliptical galaxy detected in the field with 10710 counts. ROSAT 
HRI, PSPC and Chandra data on the field were studied by Paolillo et al. (2002), in which they find that the NGC 
1404 halo is well represented by a single symmetric /3-model and follows the stellar light profile within the inner 8kpc. 
While it decreases the effective sky area available in this ObsID for serendipitous sources, our analyses do not include 
this object. 

CXOMP J033841.4-353133; Our HYDRA-CTIO spectroscopy of this bright (1710 count, r'=20, HR = -0.675) X- 
ray source shown in Figure^] reveals a BLAGN at z = 0.3597. Since it is in the background of the Fornax cluster, the 
X-ray spectrum of this quasar provides a direct measure of all the gas both neutral and ionized, along one line of sight 
through the cluster, the halos of both of its giant member galaxies, and the immediate environment of the quasar. We 
use the ACISABS model^® within CIAO Sherpa to account for the time-dependent quantum efficiency degradation of 
ACIS and fit an absorbed power-law model to the X-ray spectrum of the quasar. This yields a standard AGN slope of 
r — 1.81 ± 0.06 (90% confidence limits), with a marginal excess absorbing column of 1.6 ± 1.5 o ver the local Galactic 
column (which is 1.5, all columns in units of 10^" cm~^). Our result is consistent with those of iWaueh et al.l l|2002fl 
who find a considerable deficit of HI- rich galaxies in the centre of the cluster, but our fit further suggests that little 
ionized gas may be present, which is worthy of further investigation. 

CXOMP J033942.9-352410 is a previously known r* = 19.0, z = 1.04295 quasar FCSS J033942.9-352410. From our 
HYDRA-CTIO spectra, we measure z = 1.0417 

FCSS J033943. 1-352200 (CXOMP J033943.0-352159) is listed as an r = 17.3 NELG at z = 0.06208 

FCSS J034015.4-352850 (CXOMP J034015.4-352849) is a r* = 18.8 quasar listed at z = 0.85504, which assumes the 
bright emission line at 5225A to be Mgll. Instead, we classify this as a z = 1.74 BLAGN, because on the blue side of 
this emission line we detect the characteristic shoulder of AlIII superposed on CIII]A1909. 

Q2345+007 

The target of this Chandra image (Green et al. 2002) is a wide (7". 3) separation quasar pair Q2345-t-007A,B 
(z = 2.15), whose components have strikingly similar optical spectra (Steidel & Sargent 1991). A deep search for a 



http: / / www.astro.psu.edu / users / chartas /xcontdir / xcont.html 



38 



Green et al. 



lensing cluster revealed nothing, and differences in the absorption X-ray hardness led the authors to posit that the 
system is a binary and not a lens. Evidence for or against the lens hypothesis has been sought in optical photometric 
variability studies, allowing a search for variable objects other than the quasars. 

Three of the X-ray sources fainter than r' = 20 in this field are flagged as highly variable in Kochanski et al. and 
as such are likely to be AGN. CXOMP J234811. 6+005700: This Bj = 22 highly variable object (G5158) has 269 
counts in our ACIS image. CXOMP J234811. 9+010301 (G7676) with just 14 ACIS counts is another such source. 
CXOMP J234812. 8+005750 (G5740) has 203 ACIS counts. 

CXOMP J234813.8+005640: We classify this r* = 20.46 object as a z = 1.0416 quasar^^ based on a single broad 
emission line at 5715A. 

CXOMP J234840. 1+010753: k z = 0.71868 quasar published in the SDSS Early Data Release (2001) as 
SDSS J234840.05+01075 quasar. We include the SDSS redshift in our tabulations. 

CLJ0542.8-U00 

The intended target of this field (Harald Ebeling, PI) is a high redshift cluster. 

CXOMP J054320.2-410156.1: The extremely X-ray bright cataclysmic variable TX Col is imaged in this field. TX Col 
is a highly variable member of the in termediate polar subclass of CVs. A complete study of the Chandra observation 
of TX Col is in preparation (jSalinas , Schlcgel. & Bcckcrman 2003). 

CXOMP J054328. 5-410119: Our imaging reveals an i* — 21 object with complicated morphology, possibly 3 objects. 
A spectrum of the object shows an NELG at z = 0.3987. 

MS2137.3-2353 

This z — 0.313 cluster was chosen as a target (PI Mike Wise) for study of cooling flows. None of the X-ray sources 
in this field yielded matches in NED. 

CXOMP J213945.0-234655 is an X-ray selected z = 4.930 QSO discovered by ijSilverman et alJ I2002D from the 
ChaMP. 

CXOMP J214018. 0-234920 has an extremely blue optical spectrum, and only 1 broad line at 6735A, which wc take 
to be rest-frame MgIIA2800 at z = 1.402 ± 0.002. 



Confusingly, an object at these coordinates is listed IVeron-Cettv fc Veronir2noni^ as a ^ = 2.653 quasar with the name 2345+007C, 
but these coordinates are not near the target objects that comprise 2345+007, which in any case has z = 2.15. 



